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Abstract Fitting the spectral energy distributions (SEDs) 
of galaxies is an almost universally used technique that has 
matured significantly in the last decade. Model predictions 
and fitting procedures have improved significantly over this 
time, attempting to keep up with the vastly increased vol- 
ume and quality of available data. We review here the field 
of SED fitting, describing the modelling of ultraviolet to in- 
frared galaxy SEDs, the creation of multiwavelength data 
sets, and the methods used to fit model SEDs to observed 
galaxy data sets. We touch upon the achievements and chal- 
lenges in the major ingredients of SED fitting, with a special 
emphasis on describing the interplay between the quality of 
the available data, the quality of the available models, and 
the best fitting technique to use in order to obtain a realistic 
measurement as well as realistic uncertainties. We conclude 
that SED fitting can be used effectively to derive a range 
of physical properties of galaxies, such as redshift, stellar 
masses, star formation rates, dust masses, and metallicities, 
with care taken not to over-interpret the available data. Yet 
there still exist many issues such as estimating the age of 
the oldest stars in a galaxy, finer details of dust properties 
and dust-star geometry, and the influences of poorly under- 
stood, luminous stellar types and phases. The challenge for 
the coming years will be to improve both the models and 
the observational data sets to resolve these uncertainties. 
The present review will be made available on an interactive, 
moderated web page (sedfitting.org), where the community 
can access and change the text. The intention is to expand 
the text and keep it up to date over the coming years. 

Keywords methods: data analysis, radiation mechanisms: 
general, techniques: photometric, techniques: spectro- 
scopic, galaxies: ISM, galaxies: stellar content 

1 Introduction 

Integrated spectral energy distributions (SEDs) are our pri- 
mary source of information about the properties of unre- 
solved galaxies. Indeed, the different physical processes oc- 
curring in galaxies all leave their imprint on the global and 
detailed shape of the spectrum, each dominating at differ- 
ent wavelengths. Detailed analysis of the SED of a galaxy 
should therefore, in principle, allow us to fully understand 
the properties of that galaxy. SED fitting is thus the attempt 
to analyze a galaxy SED and to derive one or several phys- 
ical properties simultaneously from fitting models to an ob- 
served SED. This is in contrast to searching a single feature 
that could constrain a single parameter (a prominent exam- 
ple would be the use of the Ha line to derive the star forma- 
tion rate, SFR). 

The aim of this review is to present the state of the art in 
the area of fitting the ultraviolet (UV) to far infrared (FIR) 



SEDs of galaxies. It grew from and presents a summary of 
a workshop held at the Lorentz center in Leiden in Novem- 
ber 2008. As the field is extremely large we attempt here a 
somewhat novel approach to the process of writing a review. 
We consider this text as a starting point. The text will then 
be made available at www.sedfitting.org and we invite the 
community to sent us corrections and additions. We partic- 
ularly hope that this will help in covering the work that we 
might have missed. We also made no attempt (yet) to cover 
the new results of the last year, in particular from Herschel. 

Progress in many areas that affect SED fitting has been 
made recently. A major development in the last decade has 
been the advent of new observing facilities and large surveys 
at all wavelengths of the spectrum, enabling astronomers for 
the first time to observe the full SEDs of galaxies at wave- 
lengths from the X-rays to the radio. The same surveys have 
also pushed the distance of the farthest galaxy whose SED 
is amenable to study to redshifts higher than 6. At the same 
time, tools and models have been created that aim to ex- 
tract the complex information imprinted in the SEDs. Also, 
useful semi-analytic models of galaxy formation have ap- 
peared that provide realistic predictions for the properties of 
galaxies as they would be formed in the current cosmologi- 
cal standard model. Not only are astronomers exploiting the 
available data, but the next generation of surveys is now in 
the planning phase. 

For this review we concentrate on observations from the 
ultraviolet to the far infrared, including both multi-band 
photometric and spectral data. We thus treat the light emit- 
ted by stars, either directly or processed by the gas and dust 
of the surrounding interstellar medium. At wavelengths out- 
side the regime considered here, such as the X-ray and radio 
wavelengths, non-stellar processes (or at least those not di- 
rectly associated with stellar light) such as shocks, accretion 
onto compact objects, etc. dominate. While many of these 
can be associated with the star formation history of a galaxy 
(e.g. supernova rate and recent star formation), these pro- 
cesses require a higher order of complexity generally not 
considered in most current models of galaxy spectra, and 
hence we do not discuss these wavelengths in the rest of this 
review. We also do not treat the contribution of active galac- 
tic nuclei to the SEDs of galaxies. 

Our initial aim of both the workshop and this review was 
to set up a basic framework to answer the main question re- 
lating to SED fitting: considering the difficulties with the 
models, considering the limitations of the data and consid- 
ering the fitting technique, what is the true uncertainty and 
limitations on the properties that can be determined from 
galaxy SED fitting? 

This review is structured as follows: in Section [2] we re- 
view the basics of galaxy ultraviolet to infrared SED mod- 
elling, from galaxy formation to the production of and pro- 
cessing of the radiation from these galaxies. We especially 
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mention some of the current issues and main uncertainties 
of the modelling of galaxy SEDs. In Section [3] we pro- 
vide a short overview of the intricacies of assembling multi- 
wavelength SEDs. In Section [4] we present techniques and 
algorithms for SED fitting, and - most importantly - efforts 
at validating the results from the SED fitting procedure with 
independent data. Section[5]presents a review of photomet- 
ric redshift determinations, a special case of a physical prop- 
erty derived from SED fitting, as it can be compared to and 
calibrated on independently determined data, spectroscopic 
redshifts. Finally, Section [6] showcases some recent results 
from application of the SED fitting procedures, where we 
hope to emphasize the variety of problems to which SED 
fitting can significantly contribute. 



2 Modelling galaxy SEDs 

Galaxies emit across the electromagnetic spectrum. Exclud- 
ing those galaxies dominated by an accreting supermassive 
black hole at their nucleus (AGN), the ultraviolet to infrared 
spectra of all galaxies arises from stellar light, either directly 
or reprocessed by the gas and dust of the surrounding inter- 
stellar medium (ISM). Thus the UV-to-IR spectral energy 
distribution or SED contains a large amount of information 
about the stars of a galaxy, such as the stellar mass to light 
ratio, and the surrounding ISM, such as the total dust mass. 
However, to extract such information, models are necessary 
in order to connect physical properties of the galaxy with the 
observed SED. In this section we discuss such models, be- 
ginning with the stellar spectrophotometric models, moving 
on to the transfer of the radiation of these stars in a galaxy 
through the ISM, and finally how to connect these with the 
larger picture of galaxy formation and evolution. We use the 
following abbreviations for designing wavelength ranges, 
though the exact boundaries between wavelength regimes 
are not sharp: ultraviolet (UV) for /l<3500 A, optical for 
3500<^<8000 A, near infrared (NIR) 0.8</l<3 /urn, mid- 
infrared (MIR) 3<i<25 fim, far-infrared (FIR) 25<A<250 
fan, sub-mm 0.25</l<l mm, and radio A>1 mm. 

2.1 Stars 

In its simplest sense, a galaxy is a population of stars rang- 
ing from numerous, low-luminosity, low-mass stars, to the 
bright, short-lived, massive OB stars. On closer examina- 
tion, these stars are distributed in both metallicity content 
and age ranging from when the galaxy first formed to those 
newly born. The method of creating a galactic spectrum 
through the sum of the spectra of its stars is called stellar 
population synthesis and was pioneered in works by Tinsley 
( [T972) , [Searle et al.| ( fT973) and |Larson and Tinsley| ( [T978) . 
A simplification for the modelling of galactic SEDs is that 



the emitted light can be represented through a sum of spec- 
tra of simple stellar populations (SSPs) with different age 
and element abundances. Here a SSP is an idealized single- 
age, single-abundance ensemble of stars whose distribution 
in mass depends on both the initial distribution and the as- 
sumed age of the ensemble. 

There are two main methods used by current stellar spec- 
trophotometric models to compute the SEDs of SSPs: The 
first is called 'isochrone synthesis'. It uses the locus of stars 
with the same age, called an isochrone, in the Hertzsprung- 
Russel diagram and then integrates the spectra of all stars 
along the isochrone to compute the total flux. This method 
was established by Chiosi et al. ( 1988 ); Maeder and Meynet 



Chariot and Bruzual 



(1991 



and is 



( 1988) and in particular 
currently used by the majority of stellar population mod- 
els. The second uses the 'fuel consumption' approach. One 
of the problems of the isochrone synthesis method was that 
isochrones are calculated in discrete steps in time and there- 
fore phases where stellar evolution is more rapid than theses 
timesteps were not well represented (the most famous ex- 
ample of the last years being the thermally pulsing asymp- 
totic giant branch stars). Models using the fuel consumption 
theorem circumvent this problem by changing the integra- 
tion variable above the main sequence turnoff to the stellar 
fuel, i.e. the amount of hydrogen and helium used in nuclear 
burning. The fuel is integrated along the evolutionary track. 
The main idea is that the luminosity of the post-main se- 
quence stars, which are the most luminous, is directly linked 
to the fuel available to stars at the turnoff mass (for full de- 
tails, see e.g. |Buzzonifl989l |Maraston|1998| |2005l l. While 
these methods are fundamentally different in their integra- 
tion methods, most of the issues discussed here in terms of 
stellar evolution and stellar libraries apply to both. 

2.1.1 Simple stellar populations 

The spectrum (flux emitted per unit frequency per unit 
mass), L v , of a SSP of mass M, age f, and metallicity Z is 
given by the sum of the individual stars: 



Ly(t,Z)= f cf>(M\ z L v (M,t,Z). 
Jm 



(1) 



In practice the emitted light is dominated by the most 
massive, luminous stars. 

The stellar mass function, <p(M) u z, is computed from an 
initial mass function (IMF, </>o(M)) and stellar evolution, 
which describes when and which stars will stop contribut- 
ing to the SSP spectra because they end their lives either 
as Supernovae or as white dwarfs. The IMF describes the 
distribution in mass of a putative zero-age main sequence 
stellar population and is an input parameter of stellar pop- 
ulation synthesis models. The IMF is usually limited be- 
tween a minimum and maximum stellar mass (generally 
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M min ~ 0.05 - 1.0M G ; M max ~ 100 - 150M o ). Three em- 
pirical forms are most commonly used: a simple power-law 
model (Salpeter 1955; Massey 1998), a broken power-law 



HKroupa|200"T} , or a lognormal form (Chabrier 2001 ). How- 
ever whether these forms hold in all conditions and for all 
redshifts is still an open question (A good coverage of this 
field can be found in the "IMF@50" proceedings, E. Cor- 
|belli, F. Palla, & H. Zinnecker|2005| l. 

The true difficulty of calculating equation [TJ lies in the 
second part, determining the SED (L v ) of a star of initial 
mass, M, age, f, and metallicity, Z. This requires; 1) the 
computation of stellar evolutionary tracks that determine 
where a star of given stellar parameters (e.g. mass M, age t 
and abundance Z) lies on the Hertzsprung-Russel diagram or 
log g - T e tf diagram, to build up the stellar 'isochrone', and 
2) the computation or empirical building of a stellar library 
of L v with full coverage of log g, T e g, and Z to determine 
what the resulting spectrum of such a star is. 

The creation of stellar isochrones requires a large grid 
of evolutionary tracks, created by modelling the evolution 
of stars of a given initial mass and metal content. Over 
the past few decades, much work has gone into providing 
homogeneous sets of stellar tracks from different groups, 
e.g. Padova ( |Marigo and Girardi|20(T7] |Marigo et al.|2008 



Geneva ( Lejeune and Schaerer|2001[ ), Yale ( |Demarque et a 



|2004] l, MPA( W eiss and S chlattl 2008), BaSTI ( |Pietrinferni 



et al.||2009 1. For SSP modelling, the models generally run 
from the start of the main sequence (Zero Age Main Se- 
quence, ZAMS) to some end point of the star, such as a su- 
pernova or the asymptotic giant branch (AGB) phase. Origi- 
nally computed only for a solar metallicity composition and 
a few stellar masses, sets of homogeneous stellar evolution- 
ary tracks now exist for a wide range of initial masses (from 
~ 0.1M o to ~ 120 M G ; see e.g. |Girardi et al.|2000||Meynet 



|and M aeder 2005 ) and metallicities (~ 0.01 to ~ 4Z G ). In 
most stellar evolutionary modelling it has been assumed that 
for all stellar masses the elemental composition is the same 
for a given metallicity, however more recently the evolution- 
ary effects of elemental variations such as a-enhancement 
(e.g. Salasnich et al.|2000) ) or individual element variations 
(e.g. Dotter et al. 2007| ) have been investigated. However 



problems still remain in the field, with the different treat- 
ments by the different groups still giving distinct evolution- 
ary tracks even with the same inputs, as shown in figure [TJ 

While the evolutionary tracks lead to the generation of an 
isochrone, to determine a SSP spectrum, a library of stellar 
spectra is needed, covering the necessary parameter space 
in log T e ff, log g, metallicity etc. As with the evolutionary 
tracks, stellar libraries have improved significantly in recent 



years, with both fully theoretical (e.g. |Kurucz|1992[ West-| 
|era et al.|200"2"l|Smifh et al.|2002t|Coelho et al.|2003j|Martinsl 



et al.|2005j|Lancon et al.|2007[ ) and empirical, e.g. STELIB 



[20061 ICenarro et all[2007l ), Indo-US dValdes et al.|[2004| , 
ELODIE ( |Prugniel and Soubiran|200T]|Prugniel et al.|2007| ), 
HST/NGSL (Greg g~^t~aLl[2004l ) libraries covering much 
greater parameter spaces and increasing in both spectral and 
parameter resolution (see figure [2J. Unlike the optical, the 
UV still suffers from incomplete libraries which is specially 
important for fully exploiting data on high redshift galaxies 
(see e.g. Pellerin and Finkelstein|20 09 ). The question about 
which of empirical or theoretical libraries is preferable is 
a matter of the specific application (for a short review on 



both sets of libraries and their respective issues, see Coelho 
2009 ). The main benefit of the empirical libraries is that they 
are based on real stars and thus avoid uncertainties in stellar 
atmosphere structure or in opacities. On the other hand, due 
to the observational limits they cover a restricted parameter 



space biased towards Milky way compositions (see e.g. Ce- 



narro et aT]|2007[ ). Additionally, the determination of their 



fundamental parameters can be difficult for some types of 
stars and is itself based on stellar models. Conversely theo- 
retical libraries can cover a much larger parameter space and 
at any chosen resolution (see e.g. Martins et al.|2005[ ). There 
are still known problems in the comparison between the ob- 



served and theoretical stellar spectra (Martins and Coelho 
2007| ). Two specific examples for problems of theoretical 



models are incomplete line lists ( Kurucz 2005 ), problematic 
particularly at high spectral resolution, and the modelling 
of the IR emission (Lancon et al. 2007| ), which is partic- 
ularly difficult for stars in the luminosity classes I and II 
(A. Lancon, talk at workshop). The way forward may be a 
synthesizing approach, as suggest by Walc heret al.| ( [2009| ), 
aimed at using the strengths of both kinds of libraries. As 
with the evolutionary tracks, most libraries are limited to 
single compositions for a given metallicity. However re- 
cently this also has been changing, with stellar libraries ex- 
ploring abundance changes such as c-enhancement as well 
(e.g. |Coelho etal.|2007| ) . 

It is not much of an overstatement to say that the magic of 
stellar population evolutionary synthesis spectrophotomet- 
ric codes lies in interpolation. Indeed, to go from evolu- 



tionary tracks to isochrones (quoting Maeder and Meynet 
1988) "the interpolation between evolutionary tracks must 



be properly based on point of corresponding evolutionary 
status" (see also Prather|1976[ ), and similarly to go from an 
isochrone to SSP SED, the correct, often interpolated, spec- 
trum must be found for each mass bin. These stellar popula- 
tion synthesis (SPS) codes, using stellar evolutionary tracks 
and stellar libraries, then calculate Equation [TJ Besides in- 
terpolation, the challenge is to create the most homogeneous 
and most accurate set of input ingredients, interpolating in 
an appropriate manner when necessary in dependence on 
the coverage and strengths of these sets. Stellar popula- 
tion models predicting full spectra include; |Fioc and Rocca- 



( |Le Borgne etldTpOOl] ), MILES ( |Sanchez-Blazquez et al.| |Volmerange| ( [1997] PEGASE), |Bressan et al] ( [1998] used 
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in GRASIL) and ILeitherer et al. 
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Fig. 1 Evolutionary tracks of solar composition low mass stars 
(0.9 - 10 M ) demonstrating the differences between four different 
models (as labelled): MPA08 ( |Weiss and Schlattl|2008) , BaSTI04 
(with/without overshoot; Pietrinferni et al. 2009), and Padova08 
(Mari go and Girardi|2007)|Marigo et al.|2008J [Courtesy S. Char- 
lot]. 
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Fig. 2 Optical spectra from both theoretical and empirical stellar 
libraries (as labelled) demonstrating the improvement of spectral 
resolution over time with the associated improvement in library 
size [Courtesy S. Chariot]. 



(2007) (both Starburst99), [Vazdekis (1999), |Schulz et al. 



1999| ) and |Vazquez et al. 



( |2002l ), [Cervino et al.| p0"0"2"l ), |Robert et al.| pOOl] ), |Bruzual 



|and Chariot] ( |2003| GALAXE V, also commonly referred to 



as BC03), LeBorgne et al. (2004 PEGASE-HR), Maraston 



(2005, M05, based on fuel consumption theorem), |Langon 



[eTaLl ( [2008l ), |Molla et al.| ( |2009) . Currently less frequently 
used are fully theoretical stellar population models such as 
Gonzalez Delgado et al-1(|2005]l;|Coelho et al.|([2007]l. How- 



ever, as recently shown by [Watcher et al.| ([2009 ), a combi- 
nation of semi-empirical and fully theoretical models holds 
great promise for the future. As a result of the improve- 
ments in stellar evolutionary tracks and stellar libraries, as 
well as the codes themselves, stellar population synthesis 
models today can recreate broad-band UV to NIR SEDs and 
high-resolution spectra in the optical remarkably well. 

2.7.2 Validation of SSP predictions 

As SSP spectra form the basis of all fitting of galaxy SEDs 
and as the complexities of real galaxies may introduce de- 
generacies and further uncertainties in the resulting inter- 
pretations, it is of primary importance to validate directly 
the predictions of the stellar population synthesis models for 
SSPs (see e.g. Bruzual A. 2001 ). For the impact of uncer- 
tainties in the stellar parameters effective temperature, sur- 
face gravity, and iron abundance on the final SSPs see |Per-| 
|cival and Salans[ ( |2009l l. 

The ideal testbed for such validations are simple stellar 
populations occurring in nature, i.e. co-eval stellar popu- 
lations such as globular clusters (GCs), open clusters, and 
young star clusters. These objects have been used as such 
for some time (e.g. Renzini and Fusi Pecci|1988 Gonzalez 
Delgado and Cid Fernandes 2010), but unfortunately, the ex- 



act equivalence between star clusters and SSPs breaks down 
for a number of reasons: 

1) In star clusters, the stellar populations are affected by the 
dynamical evolution of the cluster, which leads to mass seg- 
regation, and evaporation of low-mass stars in GCs and, in 
young star clusters, the expulsion of gas early in their life 
time may lead to dissolution (open clusters) or to the loss of 
a significant number of stars (though see Anders et al.|2009 



for recent work on dealing with this in SSP models). 

2) Clusters also contain exotic stars (e.g. blue straggler stars, 
discussed in the following section) which influence the inte- 
grated light of the cluster but are not accounted for by most 
population synthesis models (though see Xin et al. 2007] for 
a discussion on how to account for these). 

3) Finally, star clusters only have a finite number of stars. 
In a less than 10 s M Q star cluster, the number of bright stars 
is so small, that stochastic fluctuations in the photometric 



properties of the cluster are common (Barbara and Bertelli 



1977[ Lanfon and Mouhcine 2000a] |Cervino and Luridiana 
20041 [20061 |Piskunov et al."p009| l 
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In the study of individual clusters, most of these problems 
might be alleviated by concentrating on the most massive 
specimens (W3, co Cen, starburst clusters), but these have 
the tendency to exhibit multiple rather than single (simple) 
stellar populations (e.g. Lee et al.JT999 j). Thus, one needs to 
study star cluster populations for comparison with SSP mod- 
els and to account for the influence of the stochastic fluctua- 
tions on the color-luminosity distribution. Looking the other 
way around, within a given error box for the observed colors, 
a complex distribution of possible ages is possible (Foues- 
neau et al., talk at workshop). While multi-wavelength ob- 
servations help, they do not completely eliminate the prob- 
lem. 

2.1.3 Current Issues with SSPs 

Even though significant improvements in the evolutionary 
tracks and stellar libraries have been made in the last decade, 
significant challenges remain, as some parts of stellar evolu- 
tion are only weakly understood, and hence poorly treated. 
The most important of these tend to be short lived but bright 
phases: massive stars, thermally pulsing asymptotic giant 
branch (TP-AGB) stars, extreme horizontal branch stars 
(EHB) and blue stragglers. 

Massive stars, due to their rapid evolution and short life- 
time, prove to be difficult to model and observe at all phases. 
Additional difficulties arise in that they tend to be buried by 
the interstellar material that they formed from for a large 
fraction of their lifetime, and experience high stellar winds 
(and hence strong mass evolution over their lifetime). Yet 
massive stars are a vital component of SSP modelling be- 
cause they are so luminous and can thus dominate a SSP 
spectrum, and because they give rise to most of the ion- 
izing flux and resulting nebular emission-line contribution. 
Previous modelling of massive star evolution paid particu- 
lar attention to the size of the convective core and stellar 
mass loss, yet recent theory has indicated the significant, if 
not dominant, role that stellar rotation has on the evolution 
of these stars (Meynet and Maeder 2005]>. Vazquez et al 
( 2007 ) show in their recent stellar population model that, as 
rotating stars tend to be bluer and more luminous than in 
earlier models, even the ionizing spectrum can be signifi- 
cantly altered. These differences have consequences when 
interpreting the SEDs of young galaxies, such as decreasing 
the determined mass or star formation rates. 

TP-AGB stars are short-lived, cool but luminous com- 
ponents of evolved stellar populations that tend to be more 
prominent at low metallicities. Due to the short lifetime of 
this stellar phase, as well as the inherent instability of the 
pulsations, such stars are difficult to model. Yet due to their 
relatively high luminosities (see figure [3} they can signifi- 
cantly alter the mass-to-light (M„/L) ratio of intermediate- 
age populations and it is thus important to properly include 



them. Previously, several different theoretical and semi- 
empirical recipes had been used in stellar population syn- 
thesis models leading to large discrepancies between SSP 



spectra (Vassili adis and W ood 1993 ; Maraston 1998). Atten- 
tion has been focused on these stars since Maraston (2005) 
raised this issue, leading to rapid progress in the modelling 



( |Marigo and Girardi 2007 1, largely reducing these differ- 
ences in the broad-band photometry. 

The emission of SSPs near 10 Myr is dominated by lu- 
minous red supergiants, showing that more problems exist 
in the NIR than only TP-AGB stars. At higher spectral res- 
olution in the NIR, comparisons of SSPs based on different 
libraries of synthetic stellar spectra and different isochrones 
show large residuals in the whole range from 1 to 2.5 //m 
and in particular for young to intermediate ages (A. Lancon, 
talk at workshop). Additionally, ages derived from NIR and 
optical spectroscopy are discrepant by factors of two in this 
age regime. Resolving these problems will surely lead to im- 
proved predictions of stellar population models at all wave- 
lengths. 

EHB stars represent the most luminous hot component in 
old stellar populations. Understanding these stars and im- 
plementing them in SSP codes is important because they 
could be mistaken for low-level star formation in more 
evolved, early-type galaxies. Unfortunately, the evolution of 
EHB stars is not fully understood. While HB morphology 
may be dependent upon metallicity, some metal-rich stel- 



lar populations show HB stars bluer than expected (Heber 
2008 ). The second-parameter problem with the morphol 



ogy of the horizontal branch (for a review see Catelan 2009 ) 
will need to be solved before significant progress can be ex- 
pected in this field. Meanwhile, the comparison between 
model SSP spectra and the data at old ages is affected by 
these uncertainties (Ocvirk 2 010) >. 

Blue straggler stars are, as the name suggests, stars that 
extend beyond the main sequence turn-off. Their origin is 
still unknown, though it is believed to be associated with 
binary star evolution, either through mass transfer or merg- 
ing ( |Tian et al.||2006| |Ferraro et al.1[2006| . As with EHB 
stars, blue stragglers can affect the interpretation of early- 
type spectra giving younger average ages. These stars point 
towards a fundamental limitation of current SSP modelling, 
in that effects of binary evolution are not included. While 
this causes little difference in most cases, in some situations 
(such as where blue stragglers may dominate) recent binary 



stellar population models may be more suitable (Ha n et al. 
|2007) . 

The reader is also referred to the series of papers |Conroy 
et al. ( 2009| >; Conroy and Gunn ( 20 1 0b[ l for a recent system- 



atic study of some of the uncertainties affecting SSP models. 

Yet, even with all these remaining issues, SSP modelling 
has advanced significantly in recent years, with simple Char- 
lot & Bruzual (2010, in preparation) exponentially-declining 
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SFH + burst models able to reproduce 1000's of optical 
SDSS spectra (of ~ 5A resolution) to a few percent (see 
Section [4A4l >. 

2.2 The ISM around the stars 



60 
O 




4.5 

log (T ef( /K) 

Fig. 3 Effective stellar temperature versus luminosity diagram 
showing stellar evolutionary tracks, with the problem/uncertain ar- 
eas marked. As a basis, the Padova (Marigo et al. 2008) tracks 
for 0.1 to 15 M are shown in black. At higher luminosities (blue 
curves) are evolutionary tracks for the rapidly evolving high mass 
(~ 100 M Q ) stars. The red curves show the new evolutionary 
tracks for the thermally pulsing asymptotic giant branch (TP-AGB) 
phase, while the blue and green rectangles show the regions domi- 
nated by extreme horizontal branch (EHB) and blue straggler stars 
respectively. [Courtesy S. Chariot] 



SSP models are state of the art for producing the spectra of 
stellar populations, yet they are not sufficient alone for re- 
producing the spectra of galaxies. Stars are the dominant 
power sources within galaxies (excluding AGN). However, 
the radiation from stars is absorbed and processed by the gas 
and dust that lies between the stars, the interstellar medium 
(ISM). This absorption must be accounted for when compar- 
ing SSP models with optical/UV observations and a treat- 
ment of the radiative transfer of the stellar light through the 
ISM and subsequent ISM emission is necessary if the full 
UV-IR SED is to be understood. 

While the gas and dust are in reality intermingled within 
the ISM, in practice they are often treated as separate com- 
ponents because their absorption properties have a different 
wavelength dependence. 

2.2.1 Interstellar gas 

Interstellar gas is predominantly treated as atomic in the 
modelling of galaxies. While molecular gas is clearly 
present in many galaxies, it has generally a low volume fill- 
ing factor, meaning that rarely contributes significantly to 
the overall opacity in a galaxy. It is only a noticeable opac- 
ity source in specific spectral features or in galaxies dom- 
inated by nuclear/heavily obscured sources, such as AGN 
and ultra-luminous IR galaxies (ULIRGs). Molecular gas 
emission in galaxies is predominantly seen at longer wave- 
lengths (NIR and longer) and is generally treated to arise 
mostly from "Photodissociation Regions", where the gas is 
heated by the diffuse interstellar radiation field of the galaxy. 
Although it provides insight into the molecular phase of the 
ISM, molecular emission is not considered to be a signifi- 
cant contributor to the overall SED of a galaxy (for further 
details see reviews by e.g. Young and Scoville 1991 Hol- 
lenb ach and Tielens|1997[ ) . 

Atomic gas however is the dominant opacity source in 
the extreme-UV (> 13.6 eV). It reprocesses this light into 
strong emission lines in the UV, optical and IR. It is thus 
especially important for young, actively star-forming galax- 
ies. Usually, it is assumed that all hydrogen ionizing photons 
(hv > 13.6 eV) are absorbed locally, within a small volume 
around the ionizing sources (approximately the Stromgren 
spher^]), which is appropriate due to the high opacity in this 



'The Stromgren sphere is denned as; V, = Q(H°)/(m^Q'b), where Q(H°) is 
the total number of ionizing photons, rm is the hydrogen number density, 
and ctg is the case B hydrogen recombination rate. 



s 



wavelength regime. This energy is then re-emitted in the 
hydrogen recombination lines, which correlate directly with 
the total number of ionizing photons (see e.g. |Storey and] 
[Hummer]! 995 ; Leit herer et al.|1999) . However to determine 
the emission from other elements, or to take account of both 
gas and dust absorption, full radiative transfer is needed us- 



grains, but mass dominated by large. This slope arises from 
both theory ( |Jones et aL||1996) > and matching observations 



ing photoionization codes such as CLOUDY ( Ferland et al. 



1998 ) or Mappings III ( |Groves et al.|2008| ). For further de- 
tails, see reviews by Ferland (2003 ) and Stasiriska (2007 ). 

As both the number density and absorption cross-section 
of dust is low relative to hydrogen in the EUV, dust is often 
ignored as an opacity source. However, as hydrogen absorp- 
tion is limited by the recombination rate, dust absorption be- 
comes relatively more important as the strength of the ioniz- 
ing radiation field increases, becoming the dominant EUV- 
opacity source when q(H°) > cub/k ~ 5 x 10 8 cm s _1 assum- 



ing typical values for the dust opacity, k (Dopita et al. 2002). or 



This value of the ionization parameter q is well above the 
average value measured for star-forming galaxies (see e.g. 



Kewley et al. 200 1) , meaning negligible EUV absorption by 
dust in typical H n regions (< 5%), but such high values may 
be reached within compact H n regions and AGN meaning 



dust will absorb a significant fraction of EUV photons ( Do- 
pita et al.|2002||Draine|2010) . 



2.2.2 Interstellar dust 

Interstellar dust has been a field of constant inquiry since it 
was first realized that an obscuring material existed between 
the stars and a large body of research exists on the composi- 
tion, shape and distribution of dust exists (see Draine|2003 



for a detailed review of the field, and some remaining ques- 
tions about dust). 

Most of our understanding of interstellar dust has come 
locally, from observations within our own Galaxy and the 
Magellanic clouds, and also through theoretical and exper- 
imental laboratory work. It is generally accepted that the 
grains can be considered to be composed of three differ- 
ent compositions; graphitic/amorphous carbon grains, amor- 
phous silicate grains, and polycyclic aromatic hydrocarbons 
(PAHs), which may or may not be an extension of the car- 
bonaceous grains. The former two were found to reproduce 
the observed extinction along different lines-of-sight within 
our galaxy ( |Mathis et al.| 1977 1, while the latter were added 
to explain unidentified emission bands in the mid-IR ( |Leger| 
and Puget 1984). Other forms of dust have been suggested, 
such as SiC ( Treffe rTand Cohen|1974 i, and ice is expected 
to form on grains in the coldest environments such as deep 
in molecular clouds, but generally only these forms are con- 
sidered in the SED modelling of galaxies. 

The size distribution of interstellar dust grains is thought 
to be power-law in nature, with a distribution N(a) oc a -3 - 5 or 
similar, with the average cross-section dominated by small 



( Mat his et al.|1977| |Draine and Lee|1984| [Weingartn er and| 
|Draine|200l] l. 

To obtain the optical data used for dust calculations in 
SED modelling, the size distribution and types are then con- 
volved with absorption/emission cross-sections and scatter- 
ing cross-sections and phase functions which are determined 



by both laboratory observations and Mie theory (see Draine 



|and Li|[2007al |Zubko et aLl[2004l and references within). 
PAHs are treated slightly differently as their composition is 
not fully understood, and their properties can change signifi- 
cantly with the charge of the grains, and thus have more em- 
pirical based treatments (Weingart ner and Draine|20"0T 



see 



e.g.). Altogether these form the dust models which are used 
most often in SED modelling, such as |Draine and Li| ( 2007a) 
" Zubko et aL] p004), that have been successfully compared 



with determined depletion patterns within the ISM and ob- 
servations in the UV, optical, and IR. These models are ei- 
ther used as an ensemble of individual grain sizes, or inte- 
grated to give the opacity data of dust as a whole. Of course 
empirically based laws and templates are also often used in 
SED modelling such as the Milky way extinction law and 
Calzetti law (see below). 

For the purposes of SED modelling and fitting, dust ab- 
sorption and scattering, and dust emission are often treated 
as distinct components. As the hottest dust is constrained by 
sublimation to < 2000K (corresponding to ~ 3 - 4pm peak 
emission), in practice only the scattering and absorption of 
light needs to be considered for modelling the optical-UV 
emission of galaxies. Conversely, as dust opacity strongly 
decreases with increasing wavelength, in the far-infrared 
(FIR) only dust emission needs to be considered. 

Attenuation by dust The effects of dust on the optical-UV 
light are often described by two parameters - the reddening 
and total obscuration. Reddening is the wavelength depen- 
dence of dust effects, including features, and takes account 
of the fact that shorter wavelength photons are more readily 
scattered and absorbed by dust. This is often parametrized 
by the color excess E(B - V) or the Balmer decrement 
Ha/H/3. The total obscuration is a measure of the total light 
absorbed or scattered out of our -line-of-sight by dust ei- 
ther bolometrically or in a single band and can be consid- 
ered the normalization of the reddening. This is generally 
parametrized as A(V). For relative measures correcting only 
for reddening is sufficient, however for absolute quantities 
the total obscuration must also be taken into account. This 
is especially important when the reddening is close to flat, 
i.e. only small visible effects by dust on the spectrum. 

For individual stars in the Milky Way, the Large and 
Small Magellanic clouds, extinction laws have been mea- 
sured (e.g. Cardelli et al. 1989}. However, when considering 
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a galaxy as a whole, it must be taken into account that stars 
reside at different optical depths, depending on whether they 
lie on the side of the galaxy facing the observer or averted 
from the observer, and that the stellar light can be scattered 
into the observer's line-of-sight as well as out of it. Addi- 
tionally, stellar populations of different age will have differ- 
ent extinction optical depths, and this extinction might have 
a different wavelength dependence. These issues lead to the 
concept of 'attenuation', where the complexity of the actual 
star-gas geometry is wrapped into a single attenuation law, 
now not applied individually to each star in the galaxy, but 
applied to the full spectrum of the galaxy. 

Using an attenuation law, the dust obscuration of stel- 
lar light is expressed through a screen approximation (see 
Equation[2]i, as if the dust was lying between us and the stel- 
lar population of the galaxy, with a wavelength-dependent 
reddening law (a A ). The total amount of attenuation then 
depends only upon the thickness of the screen (At), 



as bulges, disks, and dust lanes. The dumpiness of the ISM, 



I(A)obs = /stai-(/*)<? 



(2) 



The attenuation law was derived empirically for starburst 
galaxies by |Calzetti et al.| ( |1994| i; |Calzetti|fl997[ ) who fit the 
law with a simple polynomial as a function of I /A. They 
found a law much greyer than the extinction laws of the 
Milky Way and LMC demonstrating the effects of geome- 
try and mixing compared to simple extinction. Generally 
an simple power-law , a,\ oc A~ 01 , is able to reproduce the 
observed effective attenuation in galaxies (Char loTand Fall| 
12000) . 

However, a simple attenuation law cannot account for 
differential geometries and star formation histories within 
and between galaxies. This can be seen with the higher opti- 
cal depths observed for nebular emission lines relative to the 
underlying stellar continuum, indicating that the stars and 
gas that give rise to the lines and to the continuum see dif- 
ferent amounts of dust ( |Calzetti et al. |1 994} |CaIzetti| 1997) . 
These observations led to the improvement over a simple at- 



tenuation law in the approaches of Silva et al. ( 1998 1 and 



Chariot and Fall ( 2000 1, who created a more physical two- 



step model in which young stars which emit ionizing pho- 
tons are likely to be still surrounded by the clouds of gas 
and dust from which they formed. In this model all stars 
are attenuated by 'diffuse' dust in the same manner as equa- 
tion |2j However young (< lOMyr) stars undergo an addi- 
tional 'birth cloud' attenuation. In practice this means that 
the UV light and nebular emission lines associated with the 
short-lived massive stars are more obscured than the optical 
light dominated by the longer-lived stars, as observed in real 
galaxies. 



While the empirically calibrated Chariot and Fall (2000) 
model is an improvement over a simple attenuation law, it 
still does not take account of the differential dust and star 
geometries that are clearly visible in resolved galaxies, such 



both within the diffuse phase (see e.g. Kuchinski et al. 1998 



Witt and Gord on|2000[ ) and within the birth clouds (see e.g. 
|Popescu et al.|2000|[bopita et al.|2005| l, will also affect the 
resulting attenuation of galaxies. However the greatest diffi- 
culty that simple, empirically-based attenuation laws face is 
the anisotropic scattering of light by dust, as photons are not 
only scattered out of the line-of-sight, but can also be scat- 
tered into it. This can cause bluer integrated spectra than 
can be accounted for by simple attenuation laws, especially 
for face on galaxies (see e.g Baes and Dejonghe 2001b) |Frs^| 



|chera et al.|2"003)|Pierini et al.|20(M||lnoue et al.|2006 1. 

However, to take account of all these issues, proper radia- 
tive transfer (RT) calculations must be done, which require 
intensive computations. To limit these calculations several 
treatments exist, which can be broadly grouped into iterative 
methods and Monte Carlo methods (for a more detailed de- 
scription for several of the methods used in RT calculations, 
see Baes and Dejonghe 2001a). In the iterative approach, 
the light is broken up into emitted and scattered compo- 
nents, with the RT equation solved separately for each com- 
ponent, and the solution from the previous component be- 
ing used for the subsequent (i.e. directly emitted photons by 
stars, then photons scattered once by dust, photons scattered 
twice etc.) and these equations iterated to convergence (see 
e-g 



Kylafis and Bahcall 1987 Xilouris et al. 1998 1999 



Tuffs et al. 2004). Monte Carlo methods use a method closer 
to reality, where the paths of individual 'photons' are fol- 
lowed through their interactions (absorption and scattering) 
through the galaxy. The photons are emitted in a random 
direction from the sources, such as stars, and interact ran- 
domly with the surrounding ISM with a certain probabil- 
ity based on the mean free path length, and are followed 
through these scattering events until the photons escape or 
are absorbed. To build up an integrated SED of a galaxy, 
many photons must then be followed, though many treat- 
ments now exist to limit this number, such as only following 
photons which end up in the observer's line of sight (see 
e.g. |Witt et al.|1992) |Bianchi et al."p996] |Witt and Gordon] 



1996 for some early work on Monte Carlo RT in galax- 
ies). Both of these approaches are currently used, with the 
iterative quicker for given geometries, while Monte Carlo is 
more able to handle complex distributions of stars and dust 
(several existing codes are discussed in the following sec- 
tion). 

While obviously the most realistic approach, the limita- 
tion of the radiative transfer is that it requires complex cal- 
culations and thus it is not directly applicable to large sam- 
ple of galaxies. RT codes have been used to provide tem- 



plate libraries of attenuation for a range of galaxies ( Bruzual 



|A. et al.|1988"l|Ferrara et aT[ l999; Pie rini et al.|[2004) , and 
also analytic functions for the attenuation of the components 
of galaxies (i.e. bulge, disk, clumps etc., Tuffs et al.|2004 ), 
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to deal with this issue, yet these introduce several free pa- 
rameters which may be difficult to determine for unresolved 
galaxies for which only broad-band SED is available. It is 
for these reasons that a simple attenuation law is still the 
most commonly used way to account for the effects of dust 
on the UV-optical SED. 

One final note about the attenuation by dust is the sili- 
cate dust features that can appear in absorption at 9.7 and 
18 /im. These features require large optical depths to be 
observed, and thus are generally only seen in galaxies with 
strong nuclear sources (i.e. nuclear starburst/AGN). As this 
absorption occurs againstmodeled dust emission, it is usu- 
ally only modeled with a simple absorbing screen, other- 
wise it requires self-consistent radiative transfer (discussed 
in section |2".2.3| l. 

Emission by dust Dust emission in the FIR and sub- 
mm is most commonly modeled by a single black body 
(FpiR oc B^Tdust)) or emissivity-modified black body (oc 
^(Tdust)^ "A also called grey body), or a simple sum over 
a limited (2-3) number of these. The first form assumes 
that all dust is in thermal equilibrium at one temperature 
Tdust- The emissivity of dust grains is generally taken to be a 
power-law at these long wavelengths, with models and lab- 
oratory data suggesting indices ranging from (3 — 1 .0-2.0. 
Actually the /3 index is expected to be a function of both 



field and dust size and composition(see e.g. Guhathakurta 



grain size, composition and temperature (see e.g. Andriesse 



1974| Draine a nd Lee|1984l|Agladze et al.|1996||Mennella 



et al. 1 1998 with a nice discussion on the constraints on ft 
in the latter). When introducing more than one black body, 
one is generally limited by the number of wavelengths ob- 
served and the details of the model (see e.g.|Dunne and Eales 



2001 1. In general, two modified black-bodies are sufficient 



to model these wavelengths, encompassing the idea of warm 



and cold components of the ISM (see e.g. Popescu and Tuffs 



2005) . 



2002 1 Hippelein et al. 2003 and the review by Sauvage et al 



In the MIR range simple black bodies are not sufficient 
and more detailed modelling is necessary. This is due to 
strong dust (PAH) emission features and the stochastic heat- 
ing processes that become important for smaller dust grains. 
As the size of a dust grain decreases, the impingement of 
photons onto the dust grain surface becomes less frequent 
and more random, thus less statistically representative of 
the interstellar radiation field, allowing significant cooling 
between photon impacts (Figure 13 of |Draine|2003| l. Thus, 
rather than having a single temperature, the dust has a range 
of temperatures and is parametrized rather by the strength 
of the radiation field heating it. To model this one can 
use either Monte Carlo calculations simulating the arrival of 
photons and subsequent emission, or more simply one as- 
sumes and solves for a steady-state distribution of tempera- 
tures given the strength and shape of the impinging radiation 



land Draine| 1989j IDesert et al.|1990||Draine and Li|2007a| > 
Once this temperature distribution is known, it can be con- 
volved with black bodies modified by the dust emissivity in 
the MIR, including any features. 

Polycyclic aromatic hydrocarbons (PAHs) could be either 
called the largest molecular species or the tiniest dust - emit 



strong features in the MIR (see e.g. Smith et al. 2007 1. These 
features arise from specific bending and stretching modes of 
the large aromatic molecules (Bauschli cher et al.|2009| l. As 
PAH emission bands are so complex they are generally in- 
corporated into the models by either assuming a template 



form for the MIR emission features (see e.g. Desert et al. 
|1990| l or by modelling the physical processes in a way sim- 



ilar to the small dust grains (e.g. Weingartner and Draine 
|2001| |Draine and Li|2007a| ). On the whole, while aromatic 
molecules within galaxies are accepted to be the source of 
the MIR features, the typical shapes, sizes, and ionization- 
charges of these molecules are an active field of research. 

More realistic FIR dust emission models must take into 
account that the dust within the ISM of galaxies will exhibit 
a range of temperatures, from the hot dust around young 
stars and in outflows to the coldest dust in cold molecular 
cores, driven by the range of radiation fields and dust sizes. 
Such complex emission models calculate, for a given radia- 
tion field, the emission from each grain size and composition 
and then integrate over these for a given dust distribution to 
obtain the total dust emission. The largest grains are gen- 
erally considered to have a single temperature, as they will 
be in thermal equilibrium, leading to a simple distribution 
of temperatures dependent upon grain size and composition. 
In more accurate models, the smallest grains are considered 
to be stochastically heated and the temperature distribution 
of the individual grains is calculated (using, e.g., the treat- 
ment of Guhathaku rta and Draine|[l989| l. To finally calcu- 



late the IR emission from a galaxy, the distribution of dust 
masses over heating radiation field are also needed. Simpler 
IR emission models assume a functional form of dust mass 
over heating intensity; dMj = f(U)dU, with f(U) most of- 



ten assumed to be a power law (see e.g. 


Dale et al. 2001 


|Dale and Helou|2002| |Draine et al.|2007b 


). The most com- 



plex IR emission models use radiative transfer to calculate 
the radiation field distribution over a galaxy, where the dis- 
tribution of dust and stars are assumed (i.e. parameters of 
the model), and thus these models directly link the dust ab- 
sorption and dust emission. These are discussed in Section 
12331 

However, as the temperature distributions of the dust 
in the galactic ISM are dependent upon dust-gas geome- 
try and cannot be determined from optical-UV data alone, 
empirically-based templates are often used for representing 
the IR SED of galaxies, especially when IR data is limited 
due to sensitivity or confusion. These templates take dust 



11 



models as described above (i.e. multiple modified black bod- 
ies, or dust heated by a range of radiation fields) and match 
these to observed IR SEDs (or IR colors) of groups of galax- 
ies. These templates then tend to have galaxy-wide prop- 
erties such as IR luminosity or galaxy type as parameters, 
though intrinsic properties such as average interstellar radi- 
ation field intensity are also used. Well known examples 



of templates include those of Chary and Elbaz (2001 1, Dale 
and Helou] ( |2002| l, Lagache et al. (2004), and, more recently, 
Rieke et aT7| ( |2009| l. Though these templates tend to be lim- 
ited by the samples that define them, they provide a good 
alternative to models when no or very little information is 
available about the actual IR emission of a galaxy. 

2.2.3 Combining stellar and dust emission 

The full UV to IR SED of a theoretical galaxy can be created 
through the combination of the techniques and modelling 



discussed in the previous sections (2.1 -2.2.2 1. However, the 
different wavelength regimes need to be consistently con- 
nected. The simplest method is to take the energy absorbed 
in the optical-UV (see Equation[2| and to distribute it across 
the MIR and FIR, assuming simple emission properties for 
the dust, such as black bodies. This is the method used by 
|Devriendt efaX] ( fT999] l and |da Cunha et al.| pOOg] ). These 
authors attempt to strike the balance between the capabil- 
ity to model large datasets and the minimum sophistication 
necessary for a realistic model. 

To associate full UV-submm SEDs with their semi- 
analytic models (discussed in the following section) |De-| 
Ivriendt et al"1 ( f!999l ) created "STARDUST". This model as- 
sumes that stars and dust are homogeneously-mixed in the 
galaxy. The light from the stars, i.e. summed from SSPs, 
is then passed through an ISM with the amount of dust de- 
termined from a simple galaxy chemical evolution model. 
The dust-absorbed radiation is then re-emitted via a series 
of templates generated from the |Desert et al. ( 1990) model 
and fitted to observed IRAS points, parametrized by the total 
IR luminosity. 

da Cunha et al. ( 2008) > follow a similar idea, but improve 



upon this by using the Chariot and Fall (2000) recipe for 
the attenuation. They thus obtain naturally corresponding 
'birth cloud' and 'diffuse ISM' dust emission components 
over which the absorbed energy is distributed (see Figure 
|4}. The two emission components are both made up of a 
PAH template and variable grey body contributions, with the 
birth cloud emission consisting of shorter wavelength (hotter 
dust) emission. Such a model can simultaneously determine 
quantities such as stellar mass and dust mass of a galaxy, 
and provide quantitative uncertainties for all parameters (see 
While this method is quick, and hence suitable 



section 



4.1 



for comparison against large datasets, it is self-consistent 
across the two emission components only in terms of the 



total amount of radiation absorbed and re-emitted; physi- 
cal properties, such as the dust temperature or the shape of 
the emission within the components, are based on educated 
assumptions and are not constrained directly by the optical- 
UV absorption in the model. 

A very similar method was followed by (Noll et al.|2009| l 
with the CIGALE code, which uses either the Ma rastonl 



(2005) or PEGASE codes for the stellar populations and 
only a Calzetti attenuation law to attenuate the stellar light. 
The major differences lies in the use of existing empirically 
calibrated templates, such as from (Dale and Helou 2002 



rather than a free IR emission made up of several parameter- 
controlled components. 

To be properly self-consistent, the absorption and emis- 
sion must occur more 'simultaneously', such that the exact 
temperatures (including stochastic effects) of the dust caus- 
ing the absorbing can be directly calculated. Such models 
require radiative transfer calculations to be performed, such 
that the exact radiation field, or at least the heating intensity, 
is known at each point in the dusty ISM. This, along with 
assumptions about the stellar ages and distribution, and the 
dust distribution and properties can then give the full UV-IR 
SED of a model galaxy. 

The models of Efsta thiou et al.| ( |2000[ ) and |Siebenmor^| 
|gen and Kr iigel (2007 ) do this radiative-transfer calculation 
using the ray-tracing method for starburst galaxies, which, 
being dominated by young stars and their birth clouds, 
are well represented by simple spherical approximations. 
These models build upon a strong history of dust radia- 
tive transfer and emission modelling and star-formation re- 
gion modelling work to create simple models for the un- 
derstanding of the UV-submm SEDs of starburst galax- 



ies (Rowan-Robinson 1980; Rowan-Robinson and Crawford 



|1989| Rowan-Robinson|| 1 992| |Siebenmorgen and Kr uegel 
1 1992a [STebenmorgen et al. 1992b ; Rowan-Robinson and 
|Efstathiou|| 1993[ |Siebenmorgen 1993[ Krugel and Sieben- 
|morgen|[l994[ ). These works are based on the observation 
that young stars are both relatively more luminous and more 
obscured (thanks to the birth clouds) than older stars, and 
that in strongly star-forming galaxies these young stars will 
be the dominant IR (and significant bolometric) sources. In 



particular, Siebenmorgen and Krugel (2007), reduce the re- 
sults of complex modelling to a series of templates, based 
upon the physical properties of starbursting galaxies, such as 
the total luminosity, size and extinction of the star-forming 
regions and the contribution of the young stars to the total 
luminosity of the galaxy. 

Groves et al. (2008) (building on previous works; Do- 



pita et aL| |2005 2006a b), take this work a step further by 
self-consistently calculating the emission of a star-forming 
region, including the radiative transfer through the surround- 
ing gas and dust simultaneously. Like Efstathiou et al. 



(2000) they allow for the H n regions to evolve over time, us- 
ing empirically calibrated models. This model is well suited 
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for modelling starburst (star-formation dominated) galax- 
ies, where young stars and their 'birth clouds' dominate 
the emission, determining conditions such as star-formation 
rate and compactness of the gas and stars. Like |da Cunha| 
et al. (2008 1 and Siebenmorgen and Kriigel (2007} it pro- 
vides physical templates with as few parameters as possible. 
Yet, while it fits well the SEDs of star-formation dominated 
galaxies (see Figure [5j, this model is not suited for non- 
starbursting galaxies, where the distribution of the diffuse 
dust and stars must be accounted for. 

By assuming a simple molecular cloud-disk-bulge geom- 
etry (as shown in figure 151, the GRASIL model (Silva et al. 
T99"8[ |Granato et aT1|2000] > is able to account for the dif- 
ferential extinction suffered by the stars of different ages 
associated with each of these components in a galaxy. In 
addition, by varying the contribution of each component, 
galaxies from spirals to ellipticals can be modeled. Unfor- 
tunately, the more general geometry means that some parts 



(such as the gas-dust connection calculated in Groves et al 



2008 ) cannot be calculated, and also means more parameters 
are needed to define the model. As with the IGroves et al.l 
(2008 ) model, the more accurate dust calculations mean a 
longer calculation time, as compared with simpler models 
such as |da Cunha et al.| ( f20Q8] >. The GRASIL team is cur- 
rently working on speeding up their calculations for semi- 
analytic models (see following section) by the use of neural 
networks ( |Silva et al.|2010| >. 

The main issue with all models discussed above is that, 
while they take account of absorption (and emission) reason- 
ably well, they do not accurately take account of dust scat- 
tering, which, as discussed above, can make some galaxies 
appear bluer or redder depending upon inclination. This can 
be even more obvious in spatially resolved SEDs of galaxies, 
where light from stars which are obscured along our line of 
sight can be seen in reflection. However, as scattering is an 
inherently stochastic process, it is difficult to model simply 
in a galaxy, especially when multiple scatterings can occur. 

2004| l, following on from |Popescu et al 



Tuffs et al. 



(2000} and |Misiriotis et al] ( |2001| l, use the iterative ray- 



tracing radiative transfer method of Kyla fis and Bahcalr] 
(1987) to efficiently calculate the radiation field through- 
out model galaxies consisting of a stellar bulge, stellar and 
dusty disks and dusty clumps. Their resulting SEDs are then 
self-consistent across the UV-IR range. In addition, one of 
the strong benefits of radiative transfer is that the resulting 
SEDs can also be spatially resolved, and be compared to 
multi-wavelength studies of resolved galaxies, which they 
have done with edge-on galaxies such as NGC 891 ( |Popescu| 
etal.|2000| ) and NGC 5097 dMisiriotis et al.|2001| ). 

The other common approaches is to use the Monte-Carlo 
radiative transfer method to model the UV-IR SED of galax- 
ies. Existing Monte-Carlo codes that have been applied to 




\ / fj.m 



Fig. 4 IR emission of a simulated galaxy from the da Cunha et al. 
( 2008 1 model (black curve) demonstrating the individual contribu- 
tions from the 'birth cloud' dust (orange) and 'diffuse ISM' dust 
(green) [Courtesy E. da Cunha]. 




NGC 7714: SFR =9.5 M yr\ Z/Zq = l,logC=6.0,f(UCHII)=0.6 
F(old) = 3.3, A V (SB) = 1 .0 mag A v = 1 .0 mag. 



log X (f/m) 



Fig. 5 |Groves et al.| p008 ) model fit (blue curve) of the starburst 
galaxy NGC 7714 SED (black points and red curve mid-IR spec- 
tra), demonstrating the determination of physical galaxy properties 
such as star-formation rate (SFR) and metallicity (as labelled, see 
|Groves et al.|2008] for full description of parameters) [Courtesy M. 
Dopita]. 



galaxies include SUNRISE (Jonsson 2006 |Jonsson et al 
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|2010) , DIRTY( Gor don et al.|[200T] |Misselt et al.1|200l) , 
TRADING dBianchi et aTT] [19961 |2000l |Bianchi| |2008) , 
SKIRT HBaes et al.|2003l >, and RADISHE ( |Chakrabarti and] 



Whitney 2009). These are able to model arbitrary and com- 



plex geometries of dust and gas, including spiral arms, dust 
lanes, bulges and clumpy ISM. However unlike the ray trac- 
ing method, the radiation field within the galaxy is not di- 
rectly calculated (as only individual photons or photon pack- 
ets are followed). Thus dust heating and emission must be 
treated through approximations (discussed in detail within 
the papers listed above). One treatment is to integrate within 
set volumes (i.e. a grid) the amount of energy absorbed by 
dust, and to redistribute this energy over large equilibrium 
grains. In some cases (Bianchi et al. 2000), small stochas- 
tic grains are also considered (using template assumptions). 
This approach can suffer from stochastic noise if the number 
of photons used is not sufficient. A similar treatment is to 
convert the absorbed energy into a radiation field using the 
dust cross-sections, and thus with the radiation field known 
the methods described in the Section [2.2.2| can be used (see 
e.g. Misselt et al.|2001 1, though this still suffers from issues 
of stochastic noise. Another treatment is called the "dust 
temperature update method". Here, the temperature of the 
grains is updated with the absorption and emission of each 



photon (described in detail in |Bjorkman and Wood 2001 



|Baes et al.|2005) >. All these methods must iterate in the case 
of self-absorption of dust. A more efficient method iterates 
on the calculation of the radiation field density by using the 
previous estimates as a base and only calculating for the dif- 
ference at each iteration. This method will always converge 
as each iteration only adds a small amount of dust emission 
which will provide an even smaller amount of dust emis- 
sion. The new radiation field is converted to IR emission 
using models such as Dale and Helou ( 2002[ > or 



Li1< |2007a[ > (see e.g.[Juvela 2003| |Jonsson et al.|2010) . While 



Draine and 



Carlo codes require some representation of the ISM as in- 
put and are much more expensive computationally, espe- 
cially in the cases where the dust is optically thick to its own 
(IR) emission and many iterations may be required. Current 
models are also, due to resolution effects both within the 
RT and galaxy models, unable to calculate the absorption on 
both diffuse (kpc) and local (pc) scales, and thus currently 



use approximations or sub-resolution models (see e.g. Jons- 
son et al.|2010 1. Hence, while reproducing "real" galaxies, 
they cannot be directly used to fit observations of individual 
galaxies. 

In summary, the modelling of the transfer of stellar light 
through the ISM is well advanced, yet two significant chal- 
lenges still exist. The first is simply the computational effort 
needed to represent the radiative transfer accurately. Many 
of the above models are limited in their resolution to trace 
the ISM accurately, and thus need sub-resolution approx- 
imations to treat some of the coldest or hottest dust (e.g. 
SUNRISE uses the starburst templates o f|Groves et al.| 2008) 
. The second is our general lack of understanding of the 
dust composition in the ISM. Generally, dust is assumed 
to consist mainly of carbonaceous and silicate-like grains 
(such as olivine), in some power-law size distribution, (see 
e.g. |Mathis et al.|1977"| . This form is reasonably well con- 
strained by observations of extinction in the optical-UV and 
emission features in the IR (see |Draine |2003) >. Yet there are 
still open questions on shape (how ordered or "fluffy" are the 
grains, e.g. Zubko et al.|2004 1, on whether there are other 
kinds of dust, and on what formation and destruction pro- 
cesses lead to this power-law distribution of sizes (e.g. |Jones| 
et al. 1996 1. Conversely there are still spectral features asso- 
ciated with dust that are yet to be properly explained, such as 
the 2 175 A absorption feature, the diffuse interstellar bands 
in the optical, and the "E xtended Red Emission" band ob- 
served around 7000A (see 



Draine 



2003 



for a discussion on 



definitely more accurate in the treatment of dust, Monte these features and other remaining issues). 



young stars 
+ dense ISM 



diffuse ISM 
- free stars 



DISK 




BULGFI 



Fig. 6 Sketch of the geometry assumed within the GRASIL model 
(Figure 1 from |Granato et al.|2000) . 



2.3 Evolution of Galaxies 

Together, the SSP and ISM radiative transfer models of the 
previous sections are able to reproduce the full UV-sub-mm 
SED of galaxies with a reasonably high degree of accuracy 
(see e.g. |da Cunha et al.|2008j|Groves et al.|2008] >. Yet, by 
themselves, these models are inherently static. Only limited 
model assumptions about the past evolution of the galaxy 
can be introduced through the star formation history. In par- 
ticular, the ISM is rarely evolved along with the stars, and is 
presumed to be the same metallicity as the latest generation 
of stars in the SSP models. It is common form to assume 
that the dust in the ISM is a constant fraction of the metals 
within the gas, distributed in a form similar to that found in 
our Milky Way. 
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These assumptions are sufficient to reproduce the ob- 
served SEDs of real galaxies using empirically based pri- 
ors (e.g. Section |4.5| l, or multiple components (e.g. Sec- 
tion |4.4| and see Section|4]for full discussion). However to 
produce fully-theoretical SED models that are at least con- 
ceptually similar to local galaxies, one needs to fall back 
on galaxy evolution codes. There are three levels of these. 
At the innermost level are galactic chemical evolution codes, 
which, given some star-formation history and/or some "pris- 
tine" inter-galactic medium (IGM) infall rate, trace the evo- 
lution of the ISM metallicity, allowing for outflows, infalls, 
and pollution by stars (see reviews by Hensler 2008 Mat 



|teucci||2008| l. The more recent of these codes also evolve 
the dust along with the gas, taking into account the dif- 
ferent pollution rates of different elements, and the evolv- 
ing temperature/phases of the ISM (e.g. [Calura et al.|2 008). 
Once these codes have given the corresponding ISM evo- 
lution with the star-formation history (input or calculated), 
these can be associated with SSP and ISM codes to give 
a more self-consistent instantaneous spectrum of a galaxy 
(e.g. |Schurer et aT1[2009) |Conroy et ari|2010a[ ). Some of 
the main issues with these are the limited knowledge of the 
external gas losses and infalls, meaning that exact evolu- 
tion cannot be obtained, and the computational time needed 
to calculate this evolution and associate it with a spectrum, 
meaning that only specific sets of SFH or infall can be cal- 
culated at a time. 

The next scale above the chemical evolution models are 
models that evolve the whole galaxy. These models are 
based upon hydrodynamic and N-body codes that follow the 
evolution of the ISM and stars within a dark matter halo rep- 
resenting a galaxy (e.g. Springel 2005). These codes use 
empirically based relations to follow the detailed evolution, 
such as the formation of stars from gas, and the feedback 
from stars to the gas (see e.g. Tormen 1996; Cox et al. 2006). 
Containing both the stars (or "stellar particles") and the ISM 
(with known metallicity), these galaxy simulation/evolution 
codes are perfectly suited for linking with the Monte-Carlo 
radiative transfer codes such as SUNRISE (Jonsson 2006) 
or RADISHE ( |Chakrabarti and Whitney|2009| l which have 
been purposely built to create spectra and broad-band im- 
ages of these simulated galaxies. 

The outermost layer are the cosmological models. These 
trace the formation of structure in the Universe from the 
original perturbations in the cosmic microwave background 
to redshift zero, using N-Body codes to simulate dark matter 
and its gravitational interaction (see Dolag et al. 2008] for a 
review). While some of these models trace baryonic matter 
as well as the dark matter, most trace only the dark mat- 
ter due to the more complex interactions of baryonic mat- 
ter. Thus to trace the formation of galaxies within the form- 
ing dark matter halos semi-analytic models (S AMs) are used 



These models use the outputs from the dark matter sim- 
ulations and approximate the physics of galaxy formation 
within the dark matter halos by empirical relations (e.g. for 
gas cooling, star formation, AGN fueling, feedback). 

The SAMs return (and trace) the star formation history 
of each galaxy that is created, including the effects of merg- 
ers, as well as the gas content and metallicity of the gas (and 
stars). These results can be used in association with SSP 
models (as discussed in section |2.1| i to determine the stel- 
lar spectra of each galaxy. As little geometrical information 
is returned by the SAMs, associating the ISM effects on the 
stellar spectra is more difficult, especially so for the IR emis- 
sion. For the gas, most tend to use the associated emission 



lines added to the SSP models (see e.g. Leitherer et al. 1999 



Chariot and Longhetti 2001 ). For dust attenuation, a simple 



treatment taken by many is to determine the extinction as- 
suming a uniform mixing of the stars and gas in a galaxy, 
a fixed 'template' attenuation curve, and basing the optical 
depth on either empirical relations between galaxy luminos- 



ity (e.g. |Kauffmann et al.||1999] |De Lucia et al.||2004| >, or 
amount of dust in the galaxy (e.g. Guider doni and Rocca-| 
Volmerange|| 1987 Devriendt and Guiderdoni 2000). More 



advanced treatments include the use of the Chariot and Falll 
(2000) model (e.g. De Lucia and Blaizot|2007 1 or attenua- 
tion libraries like that of Ferrara et al. ( 1999 ), made for such 



purposes (e.g. |Bell et al.|2 003b 



For the dust emission, the situation is more challeng- 
ing. The simplest treatments assume that all of the radiation 
attenuated in the optical (by the above treatments) are re- 
emitted in the IR. This radiation is either distributed through 
modified Planck functions with empirically-calibrated tem- 
peratures (e.g. |Kaviani et al. 2003) > or empirically-based 



templates (e.g. Guiderdoni et al. 1998| Devriendt and 



Guiderdoni 2000 ). Yet such models do not take into account 
the strong geometrical dependence of dust heating or the 
strong variations in the spectral shape and they are clearly 
not self-consistent with the extinction in the optical-UV (see 
Section |Z23| ). 

For self-consistent SED models, the SAMs need to be 
coupled with radiative transfer (RT) calculations such as 
GRASIL, which has been done only for a few models (e.g. 
Granato et al. 2000 Lacey et al. 2008|l. However one of 



(e.g. |Cole et al . 2000; Kauff mannand Haehnelt|2000[ paF 
|ton et al.|2003)|De Lucia et al.|2004l|Somerville et al.|2008| l. 



the main strengths of SAMs is their computational efficiency 
and speed which allows the calculation of the physical pa- 
rameters of the many galaxies in large cosmological vol- 
umes and over large redshift intervals for many different 
implementations of the galaxy formation physics. Yet RT 
is computationally intensive, and severely slows the SAMs, 
meaning only relatively small volumes were investigated in 
the SAM-RT models. In addition some of the details nec- 
essary for the RT calculations are generally poorly modeled 
within the SAMs. Thus currently there is a choice between 
poorly representative but fast, or better modelling and slow 



15 



(see Fontanot et al.|[2009] for an overview). The currently 
most advanced models choose to compromise by using a RT- 



based library, empirically linked with the SAMs (Fontanot 
|et al.|2009) or even linked through artificial neural networks 
to account for the large and complex variations in galactic 
UV-IR SEDs (Silva et al., in prep). 

Models of galaxy SEDs thus exist of varying resolution 
and complexity, adapted to model everything from individ- 
ual galaxies or to large catalogs of galaxies on cosmolog- 
ical scales. While at each level of the SED models our 
knowledge of the important physical processes could be im- 
proved, SED modelling today is much more accurate across 
the wavelength range than it was even a decade ago. 



3 Constructing observed galaxy SEDs 

A major development in the last decade has been the advent 
of new observing facilities and large surveys at all wave- 
lengths of the spectrum, enabling astronomers for the first 
time to observe the full SEDs of galaxies from the UV to 
the FIR, from the local universe out to redshifts beyond 6. 
While databases such as those referenced in Table Q] make 
it tempting to simply go ahead and fit full galaxy SEDs, it 
is important to pause for a moment and review which are 
the difficulties associated with the construction of a single 
SED. Indeed, the "true SED" of a galaxy as defined in the 
models considers the sum of all photons emitted from in- 
side the volume defining the galaxy. To make the observed 
SED of a galaxy, however, this "true SED" is then filtered 
through the spectral response curve of the instruments and is 
redistributed spatially over the point spread function (PSF). 
Additionally, the measurement process not only adds noise, 
but also makes it necessary to join data from different instru- 
ments. 

The construction of multi-wavelength SEDS is a complex 
and rich subject and a dedicated review would be a welcome 
addition to the literature. In keeping with the scope of the 
present text we can here only give a very cursory treatment 
of the issue. As a starting point for further reading we sug- 
gest consulting some of the major multi-wavelength surveys 
and their overview articles provided in Table [T]below. 

3.1 Spectral response curve and resolution 

Be it in spectroscopy or in photometry, one identifies the 
SED as a series of wavelengths and associated fluxes. In 
both cases, this is only a simplification of the fact that the 
measurement process convolves the true SED with a spec- 
tral response curve, yielding a transmitted flux at an effective 
wavelength. In spectroscopy, the response curve is almost 
invariably assumed to be Gaussian, with a <j determined by 
the slit width and the dispersing device. Therefore, in prac- 
tice the distinction between the instrumental broadening and 



the broadening due to the intrinsic velocity dispersion of the 
astronomical object is not very sharp. To mimic the instru- 
mental broadening, one should first convolve with the ap- 
propriate Gaussian and then resample onto the spectral bins. 
Care needs to be taken that the wavelength calibration of 
both models and data are better than a tenth of a pixel over 
the full wavelength range (this is not always the case, see 
e.g. |Kolevaet al.|2008) . 

In photometry, the response curve is much broader and 
therefore needs to be represented with more care, i.e. tabu- 
lated as a response function. The response function in turn 
depends on the detector quantum efficiency, the instrument 
transmission and the filter in use. Photometric calibration 
and response characterization is a vital task (see Koornneef 



et al. 1986 Landolt 1992| for just two prominent examples). 



The signatures available for determination of the physi- 
cal properties of galaxies of course depend on wavelength 
and on the achieved resolution. For example, in the optical 
many of the strongest features of galaxies can be adequately 
resolved at a resolution of R=/l/A/l~2000, while the low- 
resolution part of the Spitzer IRS can easily resolve PAH 
features at R~100. However, spectroscopy is more expen- 
sive in terms of telescope time, making photometry very at- 
tractive for obtaining large samples. In the last decade, suc- 
cessful use of narrow-band filters have blurred the distinc- 
tion between spectroscopy and photometry, see for example 
COMBO- 17 HWolf et al.||2003) , COSMOS HScoville et all 
|2007) and NEWFIRM ( |van Dokkum et al.|2009| >. Narrow- 
band filters have even been used to directly measure emis- 
sion line equivalent widths (e.g. Kakazu et al.|2007] ). 



3.2 Spatial resolution, aperture bias and matching 

One of the main tasks when assembling a multi-wavelength 
SED, indeed any catalog that contains more than one mea- 
surement, is to control whether what is measured in each 
band is actually physically the same. Due to either the see- 
ing of the atmosphere or the diffraction of the telescope, the 
flux from a point-like source is re-distributed over the point 
spread function (PSF) of a width that typically depends on 
the time of observation and on the wavelength used. More- 
over, galaxies are intrinsically extended and their morphol- 
ogy may depend on the wavelength in which they are ob- 
served. 

One of the main problems in the process of matching is 
the size of the PSF. Typically, the PSF is narrowest at opti- 
cal wavelengths, while UV and IR PSFs are broader. This 
can lead to situations in which there is more than one optical 
counterpart to the UV, IR or sub-mm source. The agnostic 
way to deal with this is to simply exclude such objects from 
the sample, however, this may introduce a bias if the mul- 
tiple optical counterparts are actually physically associated. 
A more intricate, but also more uncertain, way is to redis- 



tribute the flux according to optical priors ( Guillaume et al. 
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2006 ). Finally, the use of all available information, spatial as 
well as spectral, seems to provide a promising way forward 
for multiwavelength datasets ( Roseboom et al.|20 09). 

Another problem related to the resolution of the telescope 
are aperture biases. A rather simple manifestation of this is 
that objects that are further away will be seen as smaller on 
the sky. Therefore, in order to construct comparable sam- 
ples at different redshifts, one has to adapt the size of the 
extracted aperture to the same physical size. A more com- 
plicated problem is the definition of the "total light" from 
an object. Indeed, the surface brightness profiles (SBP) of 
galaxies usually extend much beyond the threshold obser- 
vational surface brightness. In the case of specific objects, 
such as cD galaxies, these extended wings might contain 



a significant part of the total light from a galaxy (Oemler 



1976 Carter 1977 1. Different strategies have been devel- 
oped to avoid these biases, such as either integrating over 
a full model for the SBP or simply using specific apertures 
to integrate the light only inside some physical radius, but 
each method has its own problems. One of these is also 
that galaxies have different intrinsic morphology in differ- 
ent bands, thus complicating the application of consistent 
procedures, even when using data with comparable angular 
resolution. 

A particular concern with fibre spectroscopic surveys 
such as the SDSS is that the fibre apertures (3" in SDSS) 
only sample part of the object, with this fraction different at 
each redshift. The effect of this can be tested by comparing 
result of a fit to the photometry of the whole galaxy against 
the photometric fit corresponding to the area of the fibre only 
(e.g |G6mez et al.|2003[|B7Tnchmann et al.|2004) . 



3.3 Examples of multi-wavelength datasets 

We now describe three "real-life" examples of assembled 
multi-wavelength datasets. It would here be impossible to 
describe all existing datasets, we therefore chose just three 
that highlight different goals and methods. 

3.3.1 The Spitzer Local Volume Legacy - spatially resolved 
SEDs 

The Local Volume Legacy (LVL) is a Spitzer legacy pro- 
gram built upon a foundation of GALEX ultraviolet and 
ground-based Ha imaging of 258 galaxies within 11 Mpc. 
The goal of the LVL survey is to fill a vital niche in exist- 
ing multi-wavelength surveys of present-day galaxies with 
a statistically robust, approximately volume-complete study 
of our nearest star-forming neighbors. Although star forma- 
tion rates based on optical spectroscopy as well as GALEX 
ultraviolet and Spitzer infrared imaging have been measured 
for many thousands of galaxies, most currently available 
datasets are derived from flux-limited samples, and thus suf- 
fer from well-known biases against low-mass, low surface 



brightness systems. Multi-wavelength datasets that do in- 
clude such systems often only provide representative sam- 
ples of this galaxy population (e.g., SINGS; Kennicutt et al. 
|2003| >, and are thus not suitable for studies that both seek to 
probe the low metallicity dwarf galaxy regime and require 
datasets which are true to the statistics rendered by volume- 
limited sampling. LVL consolidates and builds upon recent 
Local Volume galaxy surveys which have acquired ground- 
based narrowband Ha ( |Kennicutt et al.|2008] >, GALEX ul- 
traviolet (Lee et al., in prep.) and HST resolved stellar 
population imaging ( |Dalcanton et al.| [2009), by collecting 
Spitzer IRAC and MIPS infrared imaging for a sample of 
258 galaxies derived from these programs. The collection of 
these observations enable a wealth of spatially-resolved and 
spatially-integrated studies probing present-day star forma- 
tion, chemical abundance, stellar structure, and dust proper- 
ties as well as galaxy evolution, particularly for metal-poor, 
low-mass galaxies which dominate the LVL sample by num- 
ber. The coupling of the infrared and ultraviolet data in this 



survey are explored in Section 6.2.3 



3.3.2 The Herschel ATLAS - unresolved SEDs 

ATLAS (Astrophysical Terahertz Large Area Survey) in- 
tends to provide the first unbiased survey of cool dust and 
obscured star formation in the local Universe. ATLAS will 
detect ~ 250,000 sources with a median redshift of ~ 1 
over ~ 550 square degrees of sky in five bands covering the 
wavelength range 1 10-500 yum. The five photometric bands 
cover the peak of the dust SED for local galaxies and, cru- 
cially, can probe the cold dust component (T < 20 K) which 
makes up ^ 90% of the dust mass in most galaxies ( Dunne 



and Eales 2001 ). Herschel can simultaneously provide an 



accurate measure of the bolometric far-infrared luminosity 
(related to the current star formation rate) and the total mass 
of dust (related to the gas mass). Figure [7] shows the im- 
provements expected to be made to the measurements of the 
dust mass function and the need for a multi-wavelength ap- 
proach to understanding the properties of galaxies. 

The ATLAS will be unbiased as it does not rely on prior 
detection in other wavebands. Nevertheless, the ATLAS 
fields contain > 10 5 redshifts and are the best studied fields 
of this size; they are the targets of surveys being carried 
out with GALEX, VST, VISTA, UKIRT and the South Pole 
Telescope and will be the natural targets of many future sur- 
veys, including ones carried out by DES, WISE, LOFAR 
and the two SKA precursor telescopes in the south. 

The greatest challenge for the Herschel ATLAS is likely 
to be the matching of the submillimeter sources to the cor- 
rect optical / infrared counterparts. At very low redshift this 
is fairly straightforward as the density on the sky of bright 
optical galaxies is low enough that associations are unlikely 
to be random. At intermediate-high redshifts the number of 
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Fig. 7 Top: Recent estimates of the dust mass function from 
targeted observations of an IRAS sam ple (triangles; |Dunne et al.| 
2000)) and an optical sample (stars; Vla hakis et al.||2005| l. The 
error bars at the bottom show the accuracy and range for the 
dust mass function that will be measured with ATLAS. Bottom: 
Optical-radio SEDs from GRASIL ( |Silva et al.|[T998l > for M 51 
(solid) — a typical spiral, and Arp 220 (dashed) — the archetypal 
ULIRG. Thick red lines show the SEDs at z = 0.05 where they 
have a similar optical/ultraviolet flux but very different infrared- 
radio SED. Optical / ultraviolet fluxes alone are therefore unreli- 
able as a measure of a galaxy's star formation rate in dusty galax- 
ies. Thin black lines show each SED at the redshift where it is 
just detected at 250 jim. This is z ~ 0.1 for M51 and z ~ 1 for 
Arp 220. The far-infrared-radio SEDs are now identical, while the 
optical fluxes are very different. Photometric redshifts from far- 
infrared/radio alone are totally unreliable without templates and 
optical / infrared fluxes to break the degeneracy between dust tem- 
perature and redshift. Horizontal dashes show the 5cr flux limits of 
SDSS/UKIDSS-LAS/WISE/ATLAS and LOFAR all-sky surveys. 
[Courtesy L. Dunne] 



potential matches increases dramatically and the large Her- 
schel beam (18-36") means that correct identification can- 
not simply be a matter of probability and distance. This will 
require a multi-wavelength SED modeling code which can 
'join together' the ultraviolet / optical / near-infrared portion 
of the spectrum with the far-infrared/submillimeter 'bump' 
in a self-consistent way. 

3.3.3 The SWIRE templates 



Polletta et al.| p007 ) published a set templates that combine 
the SEDs of galaxies with those of AGN and thus fill an 
important hole in many SED fitting template set^] The li- 
brary contains 20 templates including 1 elliptical, 7 spirals, 3 
starbursts, 6 AGNs, and 3 composite (starburst+AGN) tem- 
plates covering the wavelength range between 1000 A and 
1000 fim. The elliptical, spiral and starburst templates were 
generated with the GRASIL code ( |Silva et al.|1998[ >. Tem- 
plates of moderately luminous AGN, representing Seyfert 
1.8 and Seyfert 2 galaxies, were obtained by combining 
models, broad-band photometric data, and ISO-PHT-S spec- 
tra of a random sample of 28 Seyfert galaxies. Four ad- 
ditional AGN templates represent optically-selected QSOs 
with different values of infrared/optical flux ratios and one 
type 2 QSO. The QSOl templates are derived by combin- 
ing the SDSS quasar composite spectrum and rest-frame in- 
frared data of a sample of SDSS/SWIRE quasars divided in 
three groups, all, and the 25% brightest and 25% weakest 
measurements per rest-frame wavelength bin. The type 2 
QSO template (QS02) represents the SED of the red quasar 
FIRST JO 13435. 7 -093 102 ( |Gregg et al.||2002| . The com- 
posite (AGN+SB) templates are empirical templates that 
well reproduce the SEDs of the following sources: the heav- 
ily obscured BAL QSO Mrk 23 1 ( |Berta|2005| l, the Seyfert 2 
galaxy IRAS 19254-7245 South ( |Berta et al.|2003) , and the 
Seyfert 2 galaxy IRAS 22491-1808 ( |Berta|2005| l. 

3.3.4 Further examples 

While we presented above only three examples, the num- 
ber of public datasets (potentially) useful for SED fitting is 
truly staggering. Many science projects can be carried out 
without ever writing an observing proposal, the prime exam- 
ple being use of the SDSS database. This is possible thanks 
to the work of a countless number of individuals on survey 
planning, data reduction and quality insurance. Below we 
provide a table of some of the major datasets, for exploration 
and further reading. We will keep a similar table online on 
sedfitting.org and we hope to expand it in the future. 



2 The full library can be downloaded from: 
milano.inaf.it/~polletta/templates/swire_templates.html 



http://www.iasf- 
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Table 1 Some major multi-wavelength databases 



Acronym 



Full name 



Reference 
iDavis et al.|d2007k 



Website 



AEGIS All-wavelength Extended Groth Strip International Survey 

COMBO-17 Classifying Objects by Medium-Band Observations 

COSMOS Cosmic Evolution Survey 

GAMA Galaxy And Mass Assembly 

GOLDMINE Galaxy On Line Database Milano Network 

GOODS Great Observatories Origins Deep Survey 

LVL Local Volume Legacy 

MUSYC Multiwavelength Survey by Yale-Chile 

SDSS Sloan Digital Sky Survey 

SINGS Spitzer Infrared Nearby Galaxies Survey 

SSGSS Spitzer SDSS GALEX Spectroscopic Survey 

SWIRE Spitzer Wide-area InfraRed Extragalactic survey 

VVDS Vimos VLT Deep Survey 



1 1 ■ 

Scoville et al. 


(|2007 1 


Driver et al. i 


2009|> 


Gavazzi et al.l 


(2003: 


Giavalisco et al 


1 2004 


Dale etal. (2 


009| 


Taylor et al. ( 


2009i 


Abazajian et al 


(2009i 


Kennicutt et al 


fpooSi 


Treyer et al.HioToT 


Lonsdale et al. 


(2003i 


Le Fevre et al. 


(20051 



aegis.ucolick.org 
www.mpia.de/COMBO 
cosmos.astro.caltech.edu 
www.eso.org/~jliske/gama 
goldmine.mib.infn.it 
www.stsci.edu/science/goods 
www.ast.cam.ac.uk/research/lvls/ 
www.astro.yale.edu/MUSYC 
www.sdss.org 
sings.stsci.edu 
www.astro.columbia.edu/ssgss 
swire.ipac.caltech.edu 
cencosw.oamp.fr 
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Fig. 8 Four of the 20 templates from Poll etta et al.| (|2007 ), as an 
example of a set of empirically calibrated, multi-wavelength tem- 
plates [Courtesy M. Polletta]. 



4 Methods and validation of SED fitting 

In Section [2] we have described how models predict the 
SEDs of galaxies. Yet the ultimate goal of these models is 
to allow the reverse process, that is to derive the properties 
of galaxies from their observed SEDs through fitting proce- 
dures. Therefore we describe in this section several of the 
main methods used to fit model or observed template SEDs 
to multiwavelength observations of galaxies. Exactly which 
of these procedures should be used is contingent on which 
physical parameters are required and with what precision. 
In turn, this, and the accuracy with which the models repro- 
duce specific observations, determines what observational 
data are necessary. 

The following sections deal almost exclusively with stel- 
lar emission. This reflects a real dearth of specialized codes 
for carrying out the fitting process in the dust emission 
range (besides simple x 2 minimization over a parameter 
grid, e.g. Klaas et al. 2001). This state of affair is under- 
standable, however, in light of the rapid development and 
large uncertainties still inherent to our modelling of the dust 
emission of galaxies. 

4.1 Parametrizing SED models 

In the context of a SED model, one must distinguish be- 
tween independent, input parameters and derived parame- 
ters. Input parameters are those properties that are needed 
to define the SED model. Some of them may, first, not 
be measurable or, second, may not have a proper physical 
meaning by themselves. A good example for the first case 



is the shown in Figures 10 and 1 1 of da Cunha et al. ( 2008 ), 
where the width of the probability distribution functions of 
different input parameters indicates which ones are well con- 
strained, such as the global contribution of PAHs to the to- 
tal luminosity, and which are not constrained, such as the 
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temperature of warm dust in stellar birth clouds. Iglesias- 4.2 Spectral indices 



Paramo et al. (2007 their Section 3) state another example, 
namely that only a few of the many input parameters of the 
GRASIL code need to be varied in order to produce a li- 
brary of model spectra that cover a realistic range in obser- 
vational properties. Further, more systematic study of the 
relative importance of the input parameters to dust emission 
models would seem to be a important addition to the lit- 
erature. A good example of the second case is the "decay 
timescale", t»: in modelling the SEDs of early- and even 
late-type galaxies, the star formation histories (SFHs) are of- 
ten represented as falling exponential functions modulated 
by this timescale, i.e. SFR(f) oc e^'^'K However, in a hi- 
erarchical universe the SFH is likely to be more stochastic 
in form, modulated by discrete accretion events. Thus the 
width of the falling exponential that one determines from an 
SED fit is a good measure of the mean age of the galaxy. 
It should by no means be assumed that this represents the 
true SFH of the galaxy. The reason that falling exponentials 
provide reasonable fits in practice is simply that the spectral 
features of SSPs vary smoothly with time and there is thus 
considerable degeneracy between the mass contributions of 
SSPs of different age, effectively smoothing the SFHs of 
galaxies. This smoothing makes it also impossible in prac- 
tice to robustly disentangle the epoch of the formation of 
the most luminous stellar population in a galaxy from the 
timescale over which this star formation took place. 

Input parameters constitute a discrete set with a maxi- 
mum number of members determined by the quality of the 
data and the models. On the other hand the set of derived pa- 
rameters can be extended at will and not all derived param- 
eters can or need to be independent. For example, the com- 
monly used specific star formation rate (sSFR= S FR/M t ) 
and the birthrate parameter (b = S FR/(S FR),) are derived 
parameters that are very similar. It is clearly relevant to 
properly define and understand the derived physical prop- 
erties that one intends to measure. Let us again consider the 
SFR. While the SFR derived from the UV flux traces an av- 
erage SFR over the last 100 Myr, the SFR determined from 
Balmer emission lines measures the much shorter timescale 
of the ionizing stellar flux, < 10 Myr. These two measures, 
while correlated, thus need to be properly corrected before 
being able to compare them (see e.g. |Kennicutt]|1998 



Bo- 



quien et al.||2"007j |Kennicutt et~aT1[2009| |Lee et al.|[2009a 
and references therein). 

As an integrated model of the optical SED of a galaxy 
is based more or less explicitly on an equation similar to 
Equation [T] the independent and derived parameters should 
be defined in the same terms (see e.g. Wa lcher et al.|2008' 



Table 1). Actually writing down the definition helps avoid 
common confusions (such as those between metallicity Z 
and [Fe/H], or concerning the time scale over which the SFR 
is measured) and should be considered good practice. 



It has been said in the introduction that SED fitting can only 
yield useful results if the models are as or more precise than 
the effect on the data of the property to be measured. His- 
torically this was the case only for very limited wavelength 
ranges in the optical. The solution to this problem has then 
been to not fit the entire SED, but to define indices, i.e. mea- 
sure the equivalent widths, for certain absorption features. 

In the standard side-band method (Burst ein et al.| fl984; 
|Faber et aT 1985 1 a careful analysis leads to the definition of 
one central and two side bands (a blue and a red). The con- 
tinuum compared to which the equivalent width is measured 
is a linear interpolation between the average fluxes found 
in these two sidebands. Much work has gone into optimiz- 
ing the set of available indices, as well in terms of cover- 
age, as well as in terms of model precision (e.g. Rose |1984 



|Worthey et aT1[T994l |Trager et al.|[l998| [Trjpicco and Bel 



[1995] |Cenarro et aLl|2002l |Thomas et aLl|2004"l |Lee et al. 
2009b, and references therein). 

We are still waiting for standard indices in other wave- 
lengths than the optical, though see |Rix et alT] ((2004 ), Keel 



d2006| , [Maraston et al.| ( |2009> , & |Chavez et aLl ( |2009l >, and 
Lancon et aL] ( |2008| > & |CesettietaL] ( |2009| l for headway into 



the UV and NIR respectively. These would make the index 
approach a true multi-wavelength approach. 

Rogers et al. ( 2010a) i have attempted to improve on the 



classical sideband definition by introducing a "boosted me- 
dian" method. The main feature of this method is that for 
each side-band only the Lth percentile of the distribution of 
fluxes is used to determine the flux within the side-band. 
As this procedure automatically chooses the largest fluxes 
in the side-band, it will define the pseudo-continuum from 
those points that are least affected by "secondary" absorp- 
tion features, i.e. absorption features that are not part of 
the central band. It thus improves the robustness of the 
pseudo-continuum if the spectral resolution is high enough 
to avoid blending of all features. Figure |9]illustrates the dif- 
ference between the boosted median equivalent width and a 
standard side-band method for a measurement of Hy. The 
sample corresponds to a set of 14 elliptical galaxies in the 
Virgo cluster observed with a 2-2.4A resolution (FWHM) 
and with high S/N ( |Yamada et al.|2006| >. 



Index fitting can be considered a special case of SED fit- 
ting (for a fitting code see e.g. Graves and Schiavon|2008 1. 
It has the advantage of compressing the information avail- 
able in galaxy spectra into a set of discrete numbers. Much 
of what is considered secure knowledge in stellar element 
abundances and age of integrated stellar populations is still 



largely based on fitting indices (e.g. Trager et al. 2000 



IKauffmann et al.|2003aHThomas et al.|20051|Gallazzi et al.| 
[20081 |Graves etal.|2009| to cite only a few). 

The clear disadvantage of absorption line indices is that 
some of the information is lost. While indices have been 
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defined with great care to approach the ideal of representing 
one element species, in practice many small lines interfere 
in particular in the side-bands. Spectral fitting (Section |4~4| i 
can in principle use more information but it places much 
higher requirements on model accuracy. 

4.3 Principal Component Analysis 

Ideally we would like to represent a galaxy spectrum by a 
small set of continuous parameters that uniquely determine 
the best-fit spectrum. Principal Component Analysis (PCA) 
is one algorithm commonly used to derive an optimal set 
of linear components, diagonalising the covariance matrix 
of the data points to find the directions of greatest varia- 
tion. Its representation of data through a linear combination 
of independent (orthogonal) components, or eigenvectors is 
thus an alternative method to using a set of discrete SSP 
templates (Section 4.4 1. Since the convolution with trans- 



mission curves is a linear operation, these methods are as 
simple as solving a linear equation, even for photometric 
datasets ( |Connolly et al.|1999[|Budavari et al.|2000) . 

PCA has been successfully applied to astronomical spec- 
tral datasets, although not yet to photometric datasets which 
suffer the additional complication of observed-to-rest frame 
translation ( |Connolly et al.| 1995 b ). The main difficulty with 
PCA is that the interpretation of the empirically determined 
PC components in terms of physical properties is complex at 
best (though see |Wild et al.|2007|[2009||Rogers et al.|2007 



2010b). This is exacerbated by its sensitivity to outliers and 
hence out-of-the-box algorithms are of limited use for astro- 
nomical datasets. 

Recent work by Budavari et al. d2009b ) solves the prob- 



lem of reliable eigenspectra determination by an iterative 
procedure that is efficient to compute and robust in the sta- 
tistical sense. Figure 10 illustrates the comparison of three 
PCA methods on the blue optical region of 2000 randomly 
selected SDSS galaxy spectra. Each column contains (from 
left to right) the results of the classic PCA, classic PCA 
using iterative sigma clipping of the dataset, and the new 
robust algorithm. What is immediately striking in the ro- 
bust case is the following: (a) The very clean appearance 
of the nebular emission lines in the second eigenspectrum. 

(b) These are correlated with the weaker Balmer absorption 
(seen as broad wings on the narrow emission lines) and rise 
in the blue from the continua of O and B stars which are the 
dominant source of ionisation of the nebular emission lines. 

(c) The emission line in the 4th eigenspectrum is the only 
one in this wavelength range which is not a Balmer line, and 
with a higher ionisation state is sometimes attributed to the 
presence of an AGN. Without any prior physical knowledge, 
the robust PCA has separated out a line which is physically 
distinct from the others, and tied together HII emission lines 
with the O and B stars that excite them. The results are 
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Fig. 9 Comparison of the Hy F equivalent width measured with the 
side-band (SB) and the proposed Boosted Median method (Hy 2 o; 
BMC). The black dots mark two galaxies with a significant differ- 
ence between both methods. On the right-hand panels, the SED of 
these two galaxies is shown, with the slanted line representing the 
continuum defined by the SB method, and the gray line showing the 
BMC pseudo-continuum. The hatched black (gray) regions corre- 
spond to the equivalent width according to the side-band (BMC) 
method. The arrow and cross in the left panel shows the change 
in equivalent width caused by the difference in velocity dispersion 
between NGC 4387 (112 km/s) and NGC 4472 (303 km/s). [Cour- 
tesy I. Ferreras] 



-oil 5 ..ESRecJ 





SVD trimmed 




Espec 2 . 




vmU^*"** ^v*****^* 

Espec 4 . . 



3800 3900 4000 4100 4200 



Fig. 10 Comparison between the eigenspectra from a PCA of 
2000 randomly selected SDSS galaxies: left, using a simple SVD 
algorithm, the first eigenspectrum is dominated by a single noisy 
galaxy spectrum; center, using sigma clipping to remove outliers 
iteratively from the dataset, the eigenspectra are visibly less noisy; 
right, a robust PCA of the dataset, now the distinct patterns in 
the dataset are strongly visible, in particular the separation of nar- 
row nebular emission lines from broader photospheric absorption 
features. [Courtesy V. Wild] 
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clearly of more use for characterising the galaxy population 
than traditional PCA algorithms. The robust PCA algorithm 
provides a new, fast and easy to use method for the investi- 
gation of real astronomical datasets in a model independent 
manner. 

4.4 Spectral fitting by inversion 



As has been mentioned in Section [2. 1.1[ the stellar SED of 
a galaxy can be represented by a sum over the SEDs of indi- 
vidual SSPs with appropriate weights which reflect the SFH. 
As long as any complications, such as dust, can be neglected 
this is a linear problem, i.e. a matrix inversion. More gen- 
erally, inversion is the attempt to invert the observed galaxy 
SED onto a basis of independent components (SSPs, dust 
components) drawn from a SED model. Inversion is typ- 
ically used for spectral data and a big success of modern 
stellar population models and inversion codes is that we can 
now fit the models to data to better than 5% in the optical 



wavelength range (see Section 4.4.4 1 



4.4.1 Method 

Nearly all inversion methods start with assumptions that re- 
duce the complex physics of SEDs (see Eq[T| to a problem 
that can be written as a linear function of its parameters. 
Such assumptions are typically that one deals with a stellar 
system in which all generations of stars have the same metal- 
licity, i.e. £(f) = Z° and the same attenuation, i.e. T(t) = T°. 
The problem of solving for the star formation history of a 
stellar system is then equivalent to defining and minimizing 
the merit function 



X 



Fi-Z£, a k Si\t k ,ZP,T^ 



(3) 



over all non-negative a^. Here F; is the observed spectrum 
in each of n wavelength bins i, <x, is the standard devia- 
tion and a/t are weights attributed to each of M SSP models 
S i[tk,Z] of age f£ and metallicity Z. This merit function is 
linear in the a* and can thus be solved by standard math- 
ematical methods involving singular value decomposition 
(e.g. non-negative least squares, bounded least squares^] 
A particular advantage of the inversion method is that, be- 
sides the assumptions necessary to linearize the problem, no 
parametrization of the solution is necessary, in particular not 
of the star formation history. Codes that have been used for 
scientific analysis are now common and include MOPED 
( |Heavens et aTT1|2000l >, PLATEFIT ( |Tremonti et"aLp004) , 



VESPA (Tojeiro et al. 2007), STECKMAP (Ocvirk et al. 



3 It needs to be emphasized here that a correct solution involves a simulta- 
neous fit for the velocity dispersion of a galaxy. This is beyond the scope 
of the present review though. 



[2006l l, STARLIGHT ( |Cid Fernandes et al.||2005] l, sedfit 
HWalcher et alpX)^ , NBURSTS ( |Chilingarian et al.|2007 



ULySS (Koleva et al. 2009 | l. Different codes give very com- 
parable results (e.g. Koleva et al.|20 08). 

There is one problem that has to be addressed when as- 
sessing the unparameterized information content of galaxy 
spectra. All methods relying on singular value decompo- 
sition and derived algorithms inherit one of the features, 
which is that the method tends to search for the smallest 
number of templates it can use to fit the data. For galaxies, 
which presumably have a smooth SFH, this feature is a grave 
caveat concerning the significance of the recovered weights 
of each population. In particular, the recovery of realistic er- 
ror bars from a simple bootstrap algorithm is not possible, as 
the method will tend to always choose some templates over 
others, inside a range of random errors. One way to address 



(2006 



this issue is regularization, detailed in Ocvirk et al 
the STECKMAP program). 

Another robust exploration of this issue is provided in the 
code VESPA (described in detail in Tojeiro et al.|2007 i. As 
independent parameters (see |4. 1 1 > Tojeiro et al. chose 2 val- 
ues of metallicity and a logarithmic binning in age that can 
be varied between coarse (3 age bins) and fine (16 age bins). 
Dust attenuation and varying metallicity are explored by re- 
peating the fit with VESPA over a grid of parameters. Error 
bars are derived by creating noised representations of the in- 
put spectra and repeating the fit n times. The codes main 
feature in the present context is that it explores the SFHs in 
an iterative process that goes from coarse to fine resolution. 
It uses the method described in |Ocvirk et al.| ( |2006| l to es- 
timate at each step, how many parameters can be recovered 
for a linear problem perturbed by noise. The best fit will 
thus only use as many independent stellar populations as re- 
quired by the data. Nevertheless, one caveat remains, which 
is that the SFH inside each bin is fixed to be either a con- 
stant SFR or such that the light contribution of each age is 
more or less constant. When recovering parameters such as 
M/L and in those cases where the age bins are coarse, this 
will lead to an underestimate of the true uncertainty in this 
parameter, as compared to a truly non-parametric method. 

The main lesson learned from this careful analysis is that 
from spectra typical of the SDSS (S/N»;20, optical wave- 
length range), one can recover between 2 and 4 stellar pop- 
ulations described by an age and a metallicity. These val- 
ues can be improved by a factor of 2 when going either 
from S/N=20 to S/N=50 or when enlarging the wavelength 
range to include not only the optical, but also the UV range. 
Such data will be available in the near future from the 
XSHOOTER instrument at the VLT. It should be kept in 
mind, however, that while the formal constraints may be bet- 
ter, the accuracy of the results will then depend critically on 
the data calibration and model accuracy over a very large 
wavelength range. 
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4.4.2 Non-linear inversion codes 

For the optical wavelength range, an exception to the 
pre-condition of linearity exists in the form of the code 
STARLIGHT ( |Cid Fernandes et al.|[2005] ), which is based 
on simulated annealing and thus does not need linearity. 
Consistency between STARLIGHT and other codes based 
on linearity assumptions has been found. Another code 



avoiding linearity conditions was presented in Richards et al. 
( |20D9) . 

For decomposing MIR spectra into their stellar, PAH, 
dust continuum, and line emission constituents, Levenberg- 



Marquardt algorithms have been used (Smith et al. 2007 



Mars hall et al.|200~ see Section 6.2.2 1. This allows among 



other things to compute uncertainties on the derived param- 
eters. With the improvement of dust emission models and in 
particular the more widespread availability of spectral data, 
more development in the technicalities of fitting dust spectra 
can be expected. 

Finally, a possible extension of the inversion onto SSP 
models has been presented in Nolan et al. ( 2006 ) and would 
merit further exploration. Also Dye (2008} has presented 



a new method, which integrates inversion into a Bayesian 
framework and applies it to photometric data. 

4.4.3 Non-linear physics 

Dust attenuation is important in the optical and UV wave- 
length ranges. In the context of inversion, the importance 
of dust is most easily assessed by an iteration over a grid 
of attenuation values. More sophisticated schemes proceed 
in an iterative way, i.e. alternate linear inversion and a non- 
linear minimization scheme (Koleva et al. 2009 ), or use non- 



linear minimization schemes for all parameters (Cid Fernan 



des et al. 2005 1. Any estimate of the dust attenuation at 



optical wavelengths from the SED alone (i.e. using neither 
Balmer decrement, nor UV photometry, nor dust emission) 
is bound to be uncertain, in particular because no inversion 
code at present uses a physical model for different attenua- 
tions likely experienced by young and old stars. 

A very useful way is therefore to normalize the contin- 
uum of the observed and model spectra before the fit (e.g. 
|Wolf et al.||2007] |Walcher et aTJ|2009l ) or even simultane- 
ously with the fit ( Koleva et al.|2009 l. While this takes out 
any information related to the attenuation, it also allows the 
freedom from any uncertainties related to continuum cali- 
bration and can be thought of as fitting the equivalent widths 
of the absorption features only. An often under-appreciated 
caveat is, however, that this procedure throws away some in- 
formation. In particular (nearly) featureless continua (such 
as from very young stars) become undetectable, yet still af- 
fect the equivalent widths. This "featureless" continuum has 
to be appreciated in relation to the noise in the observations 



and, in an intermediate age spectrum can even include a sig- 
nificant contribution of old stars. One expression of this 
caveat is described in Section 4.6 as the bias occurring when 



fitting single SSP models to extended SFHs. 

Also the effects of forward scattering on the optical SED 
are not included in Equation[T]and cannot be linearized use- 
fully. In principle the scattering can be modelled through 
radiation transfer modelling and could thus be included in 
non-linear minimization schemes. However, radiation trans- 
fer modelling is too time-consuming to make this a practical 
possibility. 

4.4.4 Validation 

In an unpublished contribution to the workshop, Jarle 
Brinchmann (and collaborators) showed that the 10 6 spec- 
tra of galaxies in the Sloan Digital Sky Survey (SDSS) of- 
fer a rich set of tests of the current models and techniques 
for fitting SEDs in the optical region. For the first series 
of tests they use a set of 274613 galaxies from the SDSS 
DR7 with reliable determinations of the emission line fluxes. 



They fit them using the latest version of the Bruzual and 



Charlot| ( [2003[ > models (termed here CB08). In the optical 
the main improvement in these models is the switch from 
the Stelib library to the MILES library. The fitting proce- 
dure involves classical inversion using a non-negative least 
squares routine. Additionally they allow for ~ 150 A wide 
mismatches between models and spectrum using a smooth 
final correction. This is necessary because of uncertainties 
in the spectrophotometry of the data and the models. The 
result of the NNLS routine is, among others, the mean^f 2 of 
the spectrum, which can in turn be plotted into the classic 
Baldwin et ai] ( |l98l| BPT) diagram [NII]/16584A/Ha ver- 
sus [Oin]A5007A/HyS, shown in Figure 
fit is good, with a reduced x 1 around 1 : 
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Generally, the 
0.1. Regions can 
be identified, where the fit is, while still in this range, some- 
what worse. These are (A) galaxies with high star formation 
activity, high excitation and low metallicity, (B) strong, high 
excitation AGN, and (C) galaxies with high surface mass 
density (i.e. high mass) and weak emission lines. 

For cases (A) and (B), i.e., the low metallicity starbursts 
and the high-ionization AGN, they relate the main increase 
in x 2 to the imperfect masking of some emission lines and 
to the lack of nebular or AGN continua in the spectra. In ad- 
dition the BC03 models have some problems with hot stars 
(e.g. Wild et aT]|2007 1 - the preliminary version CB08 ap- 
pears to work better, due to a better coverage of hot stars in 
the spectral library (G. Bruzual, priv. comm.). 

For case (C) one of the largest mismatches in the optical 
is located, as expected, around the Mg features at 5100 A. 
This is the region where the enhancement of the abundances 
of the a elements is known to play an important role, a- 
enhancement is not covered by the CB08 models. Neverthe- 
less, in the optical the mismatch due to stellar populations is 
modest, of the order of 0.02 magnitudes. 
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The best method to validate a fitting procedure is to com- 
pare the result to independent measurements of the same 
physical property. This can be done for redshifts, as dis- 
cussed in section [5] but also for the spectra of nuclear clus- 
ters, as done by [Watcher et al. ( 2006 1 for nine nuclear clus- 
ters. The nuclear clusters have multi-aged stellar popula- 
tions and are as such similar to entire galaxies, and their to- 
tal stellar masses can be determined dynamically, indepen- 
dent from fitting the stellar spectra.. As shown in Walcher 



et al. ( 2006) 1, the stellar M/L ratios from inversion agree with 
those determined from independent dynamical modelling, 
albeit with a large scatter. This scatter reflects the general 
difficulty of constraining the oldest stellar population from 
SED fitting. Additionally, |Walcher et al.|(|2006]> also found 



that the age of the youngest population as determined from 
the spectral fitting is correlated to the measured emission 



line equivalent width, as shown in figure 12 (right). On the 



other hand, the mean mass-weighted age is not correlated, 
which is evidence that the fit result is not heavily biased to- 
wards the youngest population (in contrast to results from 
fitting single SSPs, see section [2~L2l >. 

Another careful work on validation of the spectral fit- 



ting procedures was the work by Koleva et al. (2008), who 
looked at Galactic clusters which had individual stellar spec- 
tra and full colour-magnitude diagrams to compare with the 
results of the fits to the integrated spectra (also compare this 



with section 2. 1.2 for the validation of the predictions of SSP 
models for integrated spectra). 
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Fig. 11 Classic BPT diagram of [NII]/Ha versus [OIII]/H/3, show- 
ing the distribution of SDSS emission line galaxies. The colours 
indicate the median log^ 2 of the fit of the CB08 model library to 
the SDSS fibre spectra with those emission line ratios. The upper 
and lower dashed lines indicate the Kewley et al. (2001 1 and KaufF- 
|mann et aL] ( |2003b") AGN dividing lines respectively. [Courtesy J. 
Brinchmann] 



4.5 Bayesian inference 



We refer here in particular to using the "Bayesian fitting" 
method, in which multi-wavelength SEDs are fit by first pre- 
computing a discrete set or library of model SEDs with vary- 
ing degrees of complexity and afterwards determining the 
model SED and/or model parameters that best fit the data. 
This approach is favoured in particular for multi-wavelength 
data, as the problem of solving for the physical parameters 
is not linear if effects such as dust attenuation, line emis- 
sion, and dust emission have to be taken into account. In 
addition, the influence of geometry and the physical condi- 
tions of the dust generally make the number of parameters 
very large, and with many possible degeneracies, thus mak- 
ing non-linear minimization schemes usually impractical. 

4.5.1 Method 

It is important to realise that the simple fact of pre- 
computing a set of galaxy models and afterwards determin- 
ing the one with the lowest x 1 is by itself a Bayesian ap- 
proach. By choosing which models to compute one intro- 
duces a prior, which, while possibly flat in terms of a given 
parameter, assumes that the data can be represented by that 
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Fig. 12 Comparison between the equivalent widths of emission 
lines and the mean age (left panel) or the age of the last burst (right 
panel) as determined from spectral fitting for a sample of nine nu- 
clear clusters from Walcher et al. (2006). The clear correlation is 
evidence that the youngest age picked up by the spectral fit does 
bear significance. 
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model and that parameter space. By using x 1 as a maximum 
likelihood estimator one finds the most probable model, or 
in bayesian terms, the probability of the data given that 
model. A mathematically much more rigorous description 



of the method is given in Appendix A of Kauffmann et al. 
(2003a) and currently builds the basis for Bayesian infer- 
ence in SED fitting. A more practically oriented description 



can be found in e.g. section 4 of Salim et al. (2007 1. In a 
nutshell, the method uses the fact that — assuming Gaus- 
sian uncertainties — the probability of the data (D) given 
the model (M) is given by P(D|M) oc e~ x 12 . The prior prob- 
ability distribution of the model parameters is often taken to 
be flat within a given parameter range (or flat in log space 
for larger parameter ranges). To determine the parameters 
for a specific galaxy dataset, the probabilities for all models 
are computed and then integrated over all parameters (model 
space) except the one to be derived, which yields a Proba- 
bility Distribution Function (PDF). The median and width 
of this distribution yield the parameter estimate and associ- 
ated uncertainties. The integration can be repeated and the 
PDF can be built for all parameters of the models as well 
as the observables themselves, e.g. with the aim to plot the 
model and observed fluxes in a fitted band for comparison 
purposes. 

An important step in Bayesian inference for SED fitting 
is the construction of the library and the prior assumptions 



(Section 4.5.2 1. It is equally important to synthesize the cor- 
rect observables from the models, taking into account re- 



sponse functions and redshifts (Section 3.1 1. 

Bayesian inference for SED fitting has several advan- 
tages that minimize potential sources of uncertainties: i) 
all available measurements contribute to the fit result, ii) 
the ^-correction is integrated, iii) non-linear effects such 
as from dust are accounted for as part of the models (see 
section 2.2.2 1, and thus in a self-consistent way, and finally 



iv) uncertainties on the derived parameters include measure- 
ment uncertainties as well as intrinsic degeneracies. The 
big caveats concerning this method are the sensitivity to 
the prior distribution of parameters (the "library") and, con- 
nected to this, the dependence on realistic input physics in 
the modelling. 

4.5.2 Libraries and priors 

The computation of the library is one of the essential steps 
when using Bayesian inference. It encodes our prior knowl- 
edge about the galaxies in the sample, but also our assump- 
tions about which of the physical effects we can safely ne- 
glect. Concerning photometric SEDs, a prominent case of 
neglection is the internal chemical evolution of the model 
galaxies. This is justified because even the overall metal- 
licity is not measurable from broad-band photometry alone 



(e.g. Walcher et al.|2008] l. Also the effects of forward scat- 
tering on the optical SED are often not included for compu- 
tational reasons, but could potentially be important, in par- 
ticular for starbursts ( |Jonsson|2 006). 

One assumption that is generally used is that SFRs tend 
to decrease monotonically with cosmic time in a given ob- 
ject (with the exception of star formation bursts). This as- 
sumption is particularly appropriate for early-type galaxies, 
which has led to the wide-spread use of so-called T-models, 
i.e. exponentially falling functions representing the SFH. 
However, it is not to be expected that galaxies in reality 
follow the falling exponential model. The SFHs of dwarf 
galaxies are expected to be dominated by intermittent bursts 



(see e.g. Gerola et al. 1980 1. But even for large galaxies, 



the SFH is not expected to be as smooth as suggested by the 
simple r-models. Both semi-analytic (Quillen and Bland- 



Hawthorn|2008[l and hydrodynamic modelling suggest some 



stochasticity, as discussed in section 4.1 However, due to 



the smooth variation of SED properties with age and due to 
the ensuing degeneracies, the r-models provide reasonable 
fits to observed SEDs and cannot be easily refuted. Also, 
they provide a convenient, simple parametrization of SFHs 
which explains their widespread use. 

The danger of using only r-models lies in the neglect of 
degeneracies. Let us assume, the^- 2 of a specific r-model is 
the same as that of a different model, which is the superposi- 
tion of a smooth SFH and a burst. Both models will have dif- 
ferent mean ages and M/L ratios. If they are both included in 
the library, they will thus lead to a larger error bar. However, 
if the burst model is omitted, two situations are possible: 1) 
If we indeed know that strong secondary bursts are unlikely 
(as e.g. for a sample of morphologically selected early-type 
galaxies) one will obtain a better estimate of the actual error 
bar through use of a better prior. 2) If on the other hand one 
is fitting the entire population of galaxies in a given volume, 
it is unlikely that no galaxy shows a secondary burst and 
one will thus — sometimes severely — underestimate the 
true uncertainty of the recovered parameters. Often "new", 
supposedly more precise estimates of galaxy parameters are 
published that are based on too restrictive priors. In the fol- 
lowing we showcase two examples of how to construct a 
library of star formation histories. 

A library for Gaia One of the first cases where Bayesian 
inference is applied to spectra will be presented in 
Tsalmantza, et al., (in prep). This is also a good example 
of how to create a library with a suitable prior for a specific 
application. The Gaia mission will produce spectra of 10 6 
galaxies at low resolution of R between 50 and 200. Au- 
tomatic classification of a large number of low-resolution 
spectra is therefore crucial. Classification is done by com- 
parison to a custom-created library of model galaxy spec- 
tra that simulate the properties of the Gaia spectra. This 
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library is used for training and testing of the classification 
and parametrization algorithms. 

Two libraries have been created: The first library is an en- 
tirely synthetic one. The library creation starts from the syn- 
thetic spectra of 8 typical galaxies produced by PEGASE.2 
( |Fioc and Rocc a-Volmerange 1999). The SFHs of the li- 
brary spectra are built from analytic prescriptions for either 
the dependence of SFR on the gas mass (SFR = M p gas jp1, 
late types) or for the SFH directly (SFR=£>2.e~ f/,,1 /.pL early 
types) as well as parameters for gas infall and galactic 
winds. In order to create a smooth grid, the input parame- 
ters (i.e. p\,p1, gas infall and galactic winds) were allowed 
to vary over a range of values(see Tsalmantza et al.||2007 



The crucial simplification inherent in Bell and de Jong 



2009). At this point, while of course guided by prior knowl- 
edge of the properties of galaxies, the parameter distribution 
is chosen nearly ad-hoc. 

A second library was created based on spectra from the 
SDSS (DR5, |Adelman-McCarthy and et al.|[2007| . The 
SDSS spectra were shifted to z=0 and degraded to the res- 
olution of the PEGASE.2 output spectra. They were then 
classified by comparison to the pre-existing library of purely 
synthetic spectra through simple ^-minimization. Using 
the best fit model spectrum the SDSS library spectra could 
be extended to cover the Gaia wavelength range and val- 
ues for the most significant parameters could be assigned to 
them. Mainly, however, the results of the x 2 fitting were 
used to judge how the synthetic library compares to the em- 
pirical galaxy population in the SDSS. The distribution of 
input parameters governing the SFHs obtained from the ob- 
served SDSS sample can then be used to determine a suit- 
able prior distribution for the synthetic Gaia library. This 
application is a good example of a case where the prior can 
be well established based on other observations. 

Simplifying M/L determinations A good example of how 
the appropriate choice of prior can help with constraining 
the fit is provided by the work of Bell and de Jong ( 200 1 and 
R. de Jong's talk at workshop). This work centers on finding 
the simplest way to derive the mass-to-light ratio of a stellar 
population. Let us first remember that the stellar mass itself 
is the product of a galaxies luminosity and its M/L. M/L in 
turn is a function of the star formation history of a galaxy. If 
the SFR at time t is written as *?(?) = dM(t)/dt and is defined 
over the entire Hubble time 7),, L v (t) is the luminosity of an 
SSP at frequency v, at age t, and with mass dM(t) and T v (t) 
is the mean transmission of the ISM at wavelength v and for 
the SSP with age t, then we can write 



M/L v 



jo" dm) 
jf* dtL v (t)T v (t) 



(4) 



This dependence of the M/L on the SFH means that any 
prior assumption on the SFH is also a prior assumption on 
the M/L. 



(2001 1 is thus the parametrization of the SFH. This simpli- 
fication is provided by a prior derived from a hierarchichal 
formation scenario. Using their scenario, the SFHs of galax- 
ies lead to a situation in which their intrinsic colours show 
a strong correlation with M/L. The dust attenuation vector 
is parallel to this relation, and thus the M/L ratio can be de- 
rived from a single optical colour with reasonable precision. 
Nevertheless there is a curvature in the difference between 
the M/L ratios determined from colours only and M/L ra- 
tios determined from an SED fit in the sense that for very 
blue (young) galaxies and for very old (red) galaxies, the 
colour-M/L will provide an overestimate, while for interme- 
diate colour, intermediate age galaxies the colour-M/L will 
be an underestimate. 

De Jong & Bell (talk at workshop) explore different pos- 
sible effects. When comparing a smooth SFH with a single 
burst, the M/L values have large systematic offsets, due to 
the lack of faint, old stellar populations in the single burst 
model. Two burst models do not solve the problem, in par- 
ticular if the fit is dominated by two young populations. In 
this case the two-burst model is still not complex enough, it 
again misses the faint, old population, which is crucial for 
the mass budget. The problem with the SFH template is fur- 
ther clarified by a comparison between M/L ratios as deter- 



mined in Bell et al. ( 2003a i and those derived by k-correct 
(Blan ton and Roweis 2007 1. While the scatter in the dif- 
ference is relatively small and there is no systematic trend, 
there is a systematic offset of about 0.15dex. This can traced 



to the parametrization of the SFH. Indeed, while Bell et al. 
(2003a) use pure falling exponentials, the SFH assumed in in 
k-correct peaks 6 Gyr after the beginning of star formation. 
Thus, while the resulting SED is not affected, the number 



of low-mass stars is much higher in the Bell et al. ( 2003a i 
SFHs. This uncertainty of 0.15 dex is thus inherent to any 
mass estimate from SED fitting. 

4.5.3 Validation 

In order to validate the method and its intrinsic accuracy, 
a number of studies have been performed in the last years. 
What could be termed "internal validation" is the search for 
degeneracies and systematic s in the models and the method 



itself (see e.g. |Walcher et al 2008 Wuyts et al. 2009 



Longhetti and Saracco|2 009 ) . 

Extensive testing of stellar mass determinations through 
SED fitting was done with the COSMOS dataset (Ilbe rt et al.] 
2010| l. SEDs were generated with the stellar population syn- 
thesis package of |Bruzual and Charlot| ( |2003] ) assuming a 
Chabrier ( 2003 ) initial mass function (IMF) and an exponen- 
tially declining star formation history with a timescale, t„, 
ranging from 0. 1 to 30 Gyr. The SEDs were generated on a 
grid of 51 timesteps between 0.1 and 13.5 Gyr. Dust attenu- 



ation was simulated using the Calzetti et al. ( 2000[ l law with 
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E(B - V) in the range to 0.5, with steps of 0.1. The colors 
predicted by this library of SEDs in the COSMOS dataset 
for the photometric redshift of the galaxy were compared to 
the observed ones to select the best-fitting templates. 

How the potential inaccuracy of the photometric redshifts 
affects the stellar mass derivations was investigated, and re- 
sults are displayed in Figure 13 Two spectroscopic sam- 



ples were used: the zCOSMOS bright spectroscopic sam- 
ple ( Lilly et al.|2007| l selected with a magnitude cut only at 



< 22.5 and a spectroscopic follow-up of 24 //m sources 
selected from S -COSMOS ( |Sanders et aLl[2057l ). These 
sources are typically dusty starbursts with 18 < < 25. A 
median difference smaller than 0.002 dex is found between 
the photo-z and the spectro-z stellar masses, for both spec- 
troscopic samples. 

It is interesting to note that the stellar synthesis models 
have the same effect for both selections, as shown in Fig- 
This is not the case when studying the impact of 



ure 
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the attenuation law. The MIPS selection consisting prefer- 
entially of dusty galaxies, the difference due to the choice of 
attenuation (the Cal zetti et al.|2000| law in a case, the |Char-j 
lot and Fall 2000 in the other) is amplified for this sample. 



The median difference in stellar masses is -0. 1 dex for the 
optical selection and -0.24 dex for the infrared selection. 

The impact of the choice of stellar synthesis model on 
the stellar masses determinations was also examined. Fig- 
ure 13 presents the difference between stellar masses com- 



puted using the Bruzual and Chariot (2003 ) models and the 
new models of Chariot and Bruzual (2007, private commu- 
nication) including a better treatment of the pulsing asymp- 
totic giant branch. A median difference of 0.13-0.16 dex 
and a dispersion of 0.09 dex is observed. Pozzetti et al. 



( 2007| > found a similar difference (0.14 dex) when compar- 
_ d2003l and|Maraston|((2005l) mod- 



mg 
els. 



Bruzual and Chariot 



Walcher et al. 



(2008) show that this offset changes 
with the measured sSFR of the galaxy, in the sense that the 
masses determined using BC03 for objects with an interme- 
diate sSFR (-12 < log(sSFR) < -9) are, in the mean, more 
massive by 50% than those determined using CB07. These 
are the objects where the TP-AGB stars are most likely to 
contribute significantly to the light in the NIR bands. Ob- 
jects with lower or higher sSFRs do not show this offset in 
IWalcher et al.1 ( |2008"l ). 

|Ilbert et aT (20101 concluded from this study that the use 
of photometric redshifts is only a weak source of uncertainty 
when deriving stellar masses for a sample of galaxies. The 
uncertainties are dominated by the uncertainties in the SED 
modeling itself, thus one has to be very cautious about the 
interpretations when selecting samples where a specific type 
of model is preferred. 

External validation of the method is possible for one pa- 
rameter in particular, namely the SFR. Using the entire SED 
to derive the SFR was considered less precise than other 



(photo-z - spectrp-z) 

(Calzetti - CharlotJ) 
--- (BC03 - CB07) 



(photo-z - apectr 

(BC03 - CB0T) 


Ui 

MIPS selected 


j 
...... .nl». 





-O.fl -0.4 -0.3 



-0.8 -0.4 -0.3 0.2 0.4 0.6 
A(log (M)) 



Fig. 13 Histograms of the difference between two stellar mass de- 
terminations, for the zCOSMOS bright sample selected at i + AB < 
22.5 (right) and for a sample of sources selected at 24 fim in the 
S-COSMOS MIPS sample (right). The black histograms show the 
difference in stellar masses when using the photometric or the spec- 
troscopic redshifts. The thin black dashed lines are gaussian distri- 
butions with a = 0.02 (left panel) and a = 0.03 (right panel). The 
green dashed lines show the differences between the stellar masses 
computed with Bruzual and Chariot ( 2003 1 and Chariot & Bruzual 
(2007, private communication). The magenta dotted lines show the 
differences between the stellar masses computed using the |Calzetti] 
|etal.| l |2000t and |Charlot and Fall| ( [2000) attenuation laws. The red- 
shifts were set to the spectro-z values in the two last cases. System- 
atic uncertainties due to the models dominate the errors introduced 
by the photo-z, at least in the magnitude/redshift range explored 
with our spectroscopic samples. [Courtesy H. Aussel] 
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methods (based on specific tracers) in Kennicutt (1998 ). The 
availability of large samples with UV, optical and NIR pho- 
tometry have lead to a reappraisal of the method in the last 
years. |Salim et aL (2007 ) show a detailed comparison based 
on the data from the SDSS and the Galex satellite. They de- 
rive SFRs using the multi-wavelength SED. They then com- 
pare these to the SFRs derived by Brinchmann et al. ( 2004), 
which used all emission lines that can be usefully measured 
in the SDSS spectra and detailed modelling of the emission 
line spectrum to derive dust-corrected, emission-line based 
SFRs. As shown in their figure 6, the SED-derived SFRs 
show in general a satisfying agreement. A similar test was 
performed in |Walcher et aL] ((2008]) on a VVDS-SWIRE- 
GALEX-CFHTLS sample, where the SED SFR was com- 
pared to SFRs derived from the [OII]/l3727A emission line. 
While the latter are much less accurate than the emission 
line SFRs from Brinchm ann et al.| ( |2004[ ), the agreement be- 
tween SED and emission line SFR was shown to be satisfy- 
ing out to a redshift around one (their figure 8). An advan- 
tage of obtaining SFRs from broad-band SEDs compared to 
emission-lines is the time needed to obtain sufficient S/N 
in the spectra to get the lines. Nevertheless, an interest- 
ing residual correlation of SFR difference with stellar mass 
is found in SED fitting, which remains yet to be fully un- 
derstood. Another interesting problem exists in early-type 
galaxies, in which no Ha emission can be measured, yet UV 
emission is commonly found (see e.g. Salim et al. 2007 Fig- 
ure 3). This probably is the cause why some objects with no 
measurable emission lines have relatively high SED SFRs. 

In order to test the SFR determinations from SED fitting 
at z~l, Salim et al. (2009) extended previous work to in- 
clude the mid-IR dust emission. The consensus is that the 
mid-IR emission is heated by the UV, and hence traces the 
emission of young stars (< 100 Myr) and it therefore pro- 
vides a good means to measure the SFR. Deep MIPS 24 /mi 
imaging exists for the Extended Groth Strip (EGS), as well 
as redshifts assembled in the context of the DEEP2 survey. 
The rest-frame 12 /im emission was converted to a total dust 
luminosity using the |Dale and Helou| ( |2002| ) prescriptions. 
Using the dust attenuation determined from fitting the UV- 
optical SED one can calculate the dust-corrected UV and 
B-band luminosities. These correlate well with the 24/mi 
luminosity as seen in figure 14 Turning the IR luminosity 
into a SFR using the relations of [Kennicutt ( 1998), the best 
correlation was found for an SED-SFR measured over the 
last 1-3 Gyr. That the mid-IR may be heated by older stars 
recently received support in a study by Kel son and Holden] 
( [2010] ), which predicts that C-type TP-AGB stars (-1.5 Gyr 
old) can produce a large fraction of the mid-IR flux. 

Nevertheless, for some very red objects in their sample, 
some SED SFRs seriously underpredict the SFR determined 
from the IR emission. These galaxies have a dominant old 
population, where the IR may have little to do with SF. In- 
deed, for the case of early-type galaxies Johnson et al. (in 
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Fig. 14 Comparing UV and optical luminosities to the IR lumi- 
nosity determined from the 24yum flux. The offset (A) and disper- 
sion (cr) around the mean relation of the galaxy sample are marked 
in the lower right corner. Figure from Salim et al. (2009 ) [Courtesy 
S. Salim]. 
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prep.) find that it is important to consider the fraction of 
the stellar emission absorbed redwards of 4000 A in order 
to obtain a good prediction for the measured FIR emission. 
Neglecting this emission leads to an underprediction of the 
FIR emission by factors of up to 5. 

Stellar mass is the second parameter that can in princi- 
ple be calibrated empirically, with the hope of deriving the 
normalization of the M/L ratio, in other words the choice of 
IMF. However, uncertainties related to technical questions 
concerning the SED fitting and the dynamical mass mea- 
surements are still too high to derive very accurate overall 
IMF normalizations (see e.g. Cappellari et al. 2006; S alucci| 
et al.|2008 >. Careful attention should also be paid to issues 
such as the contribution of dark matter, recycled gas, and 
other dark components. In the formulation of Bell and de 



prescription for chemical enrichment (CXP). All these mod- 



Jong (2001 ), the slope of the colour-M/L relation is indepen- 



dent of the IMF, but the normalization depends on it. For this 
reason Bell and de Jong ( 2001| l used the maximum disk ap- 
proximation to provide an upper boundary of the total mass 
in stars allowed in spiral galaxies. They thus confirm that 
the Salpeter IMF leads to too high M/L ratios and they adapt 
the normalization of the colour-M/L relation to be 0.15 dex 
below Salpeter. Thus a Kroupa or Chabrier IMF seem to be 
good choices (see section 2.1.1) . Variations between stellar 
population models and in dependence on the available data 



range also need to be considered carefully (see e.g. van der 
Wei et aL|2006 i. As many more independent mass tracers 
exist, which require careful assessment, the interested reader 
is pointed to de Jong & Bell (in prep.) for further discussion 
of this topic. 

4.6 Method-independent caveats 

A good example of how the quality of the available models 
and data influence the use of tracers and the precision with 
which physical properties can be recovered is given by re- 
cent developments in the use of index fitting. Historically, 
early-type galaxies have been analyzed using simple stellar 
populations models as templates, effectively using the prior 
assumption that early-type galaxies were created in one sin- 
gle burst of star formation. However, detailed spectroscopic 
studies ( Trage r~et al.|2000 1 as well as near-ultraviolet pho- 
tometry of early-type galaxies ( Ferreras and Silk||2000 Yi 
et al. 2005 Kaviraj et al.|2007 1 have confirmed the presence 



of hot stars in early-type galaxies. 



Figure 15 illustrates that it is now possible and necessary 



to go beyond SSPs for early-types (see also e.g. Serra and 



Trager 2007). The marginalized age distribution is shown 
for four galaxies corresponding to three different models 
(from left to right): Simple Stellar Populations (SSP), a 2- 
burst model consisting of a superposition of an old and a 
young SSP (2BST), and a composite model assuming an ex- 
ponentially decaying star formation rate, including a simple 



els combine the population synthesis models of Bruzual and 
|Charlot| ( |20"0"3| >. The figure shows that SSPs give mutually 
inconsistent age distributions, whereas composite models 
such as 2BST or CXP give a more consistent picture. Notice 
how lower mass galaxies such as NGC 4489 or NGC 4239 
(top panels) give better fits for a 2 burst scenario, whereas 
higher mass galaxies (NGC 4464, NGC 4387; bottom pan- 
els) are better fit by a smooth star formation history. 



A similar analysis is presented in Trager and Somerville 



(2009), who analyze mock data from semi-analytic mod- 
els in standard ways, in particular computing the "SSP- 
equivalent" ages and metallicities. They quantify that the 
SSP-equivalent age is poorly correlated with the mass- 
weighted or light-weighted average ages. The SSP-ages 
tend to be younger, biased to the star formation in the last 
0.1-2 Gyr. This has in particular the effect of exaggerat- 
ing the signature of downsizing. On the other hand, the 
SSP-equivalent metallicity is mostly equivalent to the light- 
weighted metallicity. This serves as an important reminder, 
that the prior assumption that one puts into the analysis are 
very important for the resulting measurement. In the case 
of using SSP-equivalent ages, the assumption, or prior, is 
that galaxies are single-age entities. This inevitably leads to 
important biases in the determined parameters. 

Among the important caveats of SED fitting, the issue of 
the age-metallicity degeneracy cannot remain unmentioned. 
This degeneracy is particularly important in the optical. In a 
nutshell, age effects on colors or line strengths can be mim- 
icked by a change in metallicity. This is especially important 
in the old populations found in early-type galaxies. Equally 
important and generally well-understood are the effects of 
using the optical SED only to determine the attenuation. Ba- 
sically this is a very difficult undertaking and under almost 
any circumstance, measurements of the Balmer decrement, 
the FIR emission or the UV slope are essential to a reliable 
determination of the attenuation. 

In principle it should be possible and informative to com- 
bine information from the photometric SED and spectral in- 
formation dGallazzi and Bell|2009l|Lamareille|200"9] l. While 
widespread use of this combination is hampered by possible 
systematic differences in the results (e.g. |Wolf et aT]|2007 



|Schombert and Rak os 2009), general agreement between 
both types of information has been claimed at least for stel- 
lar masses (Drory et al. 2004 1. Brinchmann et al. (Talk 



at workshop) showed the results of a systematic compari- 
son between photometric and spectral SED fitting. Compar- 
ing formal errors, they find no significant systematic offset. 
However, for actively star forming galaxies it might be bet- 
ter to use multiwavelength, broad-band photometry than the 
detailed spectrum information. The reason is, as mentioned 
in Section 4.4.3 that very young stars tend to "hide" in the 



normalized optical part of the spectrum, however they are 
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readily visible as blue continuum in the multi-wavelength 
SED. On the other hand, spectra are better at picking up re- 
cent bursts through their Balmer lines (e.g. |Wild et al.|2007| 
2009). Emission lines are a potential problem for fitting 



photometric SEDs. A quick back of the envelope calcula- 
tion, however, shows that they produce a maximal offset in 
r-i colour of 0.1 at an equivalent width of 100 A. Thus, ELs 
are a minor issue in the broad-band colours when a full SED 
is available and of limited concern for normal z~Q galaxie^] 
Galaxies with significant redshifts can in principle be 
treated the same way as local ones, subject to two main 
caveats though: (1) we currently do not have the same kind 
of information available for large samples of galaxies at red- 
shifts above 1 as for local galaxies. (2) Galaxies at high 
redshifts may have been significantly different from todays 
galaxies, e.g. concerning their typical star formation histo- 
ries, their metallicities or their gas content. As local analogs 
are rare or lacking, SED models are less well calibrated 
and may be subject to considerable systematic uncertain- 
ties. These need to be explored in detail, which is cur- 
rently only possible through semi-analytic models (see e.g. 
Schurer et al.|20 09 ). The design limits of SED fitting codes 
should thus be kept in mind when quoting results and - in 
particular - errorbars on high redshift properties. 



5 Results of SED Fitting: Photometric redshifts 

A special case of analyzing SEDs of extragalactic sources 
is the problem of redshift estimation, a topic that is usually 
refered to as photometric redshifts (hereafter photo-z). This 
problem is distinct from all other estimates of physical prop- 
erties because independent and more precise measurements 
of the same property are available for large samples in the 
form of spectroscopic redshifts. The method can thus be 
tested extensively and even calibrated empirically. It is also 
one of the earliest forms of SED fitting, having been sug- 
gested as a manner to go beyond the limits of early spec- 
troscopy (Bau m|1957] l. 

For a working definition, Koo ( 1999| l suggests the fol- 
lowing: "photometric redshifts are those derived from only 
images or photometry with spectral resolution A/AA < 20. 
This choice of 20 is intended to exclude redshifts derived 
from slit and slitless spectra, narrow band images, ramped- 
filter imager, Fabry-Perot images, Fourier transform spec- 
trometers, etc." This definition leaves room for a wide vari- 
ety of approaches that are actively being explored by mem- 
bers of the community. While today most studies build on a 
set of magnitudes or colors, recently other observables have 



and Gunn ( 2008} . However, all methods depend on strong 
features in the SEDs of the objects, such as the Balmer break 
or even PAH features (Negre llo et al.|2009) . 

Traditionally, photometric redshift estimation is broadly 
split into two areas: empirical methods and the template- 
fitting approach. Empirical methods use a subsample of the 
photometric survey with spectroscopically-measured red- 
shifts as a 'training set' for the redshift estimators. This 
subsample describes the redshift distribution in magnitude 
and colour space empirically and is used then to calibrate 
this relation. Template methods use libraries of either ob- 
served spectra of galaxies exterior to the survey or model 
SEDs (as described in Section [2|. As these are full spectra, 
the templates can be shifted to any redshift and then con- 
volved with the transmission curves of the filters used in the 
photometric survey to create the template set for the redshift 
estimators. 

Both methods then use these training sets as bases for 
the redshift estimating routines, which include ^-fitting and 
various machine learning algorithms (e.g. artificial neural 
networks, ANNs). The most popular combinations are y 1 - 
fitting with templates and machine learning with empirical 
models. For a review of the ideas and history of both areas, 
see |Kool ( [T999) . 

The preference of one over the other is driven by the limi- 
tations of our understanding of the sources and the available 
observations. Template models are preferred when explor- 
ing new regimes since their extrapolation is trivial, if poten- 
tially incorrect. Empirical models are preferred when large 
training sets are available and great statistical precision is 
required. Here we review these techniques and estimators, 
concentrating predominantly on the template method which 
is closer to the idea of SED fitting as discussed in the previ- 
ous section. 

5.1 Methods 

5.1.1 Empirical techniques 

Early on, the first empirical methods proved extremely pow- 



erful despite their simplicity (see e.g. Connolly et al. 1995a 



been utilized with good success, e.g., in the work by Wray 



4 But recall that equivalent widths are proportional to (1+z) and that the 
sSFR also tends to increase with redshift. 



|Brunner et al.|1997||Wang et al.|1998) . This was partly due 
to their construction, which should provide both accurate 
redshifts and realistic estimations of the redshift uncertain- 
ties. Even low-order polynomial and piecewise linear fitting 
functions do a reasonable job when tuned to reproduce the 
redshifts of galaxies (see e.g. |Connolly et al.|1995a i. These 
early methods provided superior redshift estimates in com- 
parison to template-fitting for a number of reasons. By de- 
sign the training sets are real galaxies, and thus suffer no 
uncertainties of having accurate templates. Similarly as the 
galaxies are a subsample of the survey, the method intrin- 
sically includes the effects of the filter bandpasses and flux 
calibrations 
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One of the main drawbacks of this method is that the 
redshift estimation is only accurate when the objects in the 
training set have the same observables as the sources in 
question. Thus this method becomes much more uncertain 
when used for objects at fainter magnitudes than the train- 
ing set, as this may extrapolate the empirical calibrations in 
redshift or other properties. This also means that, in prac- 
tice, every time a new catalogue is created, a corresponding 
training set needs to be compiled. 

The other, connected, limitation is that the training set 
must be large enough that the necessary space in colours, 
magnitudes, galaxy types and redshifts is well covered. This 
is so that the calibrations and corresponding uncertainties 
are well known and only limited extrapolations beyond the 
observed locus in colour-magnitude space are necessary. 

The simplest and earliest estimators were linear and poly- 
nomial fitting, where simple fits of the empirical training set 
in terms of colours and magnitudes with redshift were ob- 
tained (see e.g. |Connolly et al.|1995a i. These could then be 
matched to the full sample, giving directly the redshifts and 
their uncertainties for the galaxies. Since then further, more 
computational intensive algorithms, have been used, such as 
oblique decision tree classification, random forests, support 
vector machines, neural networks, Gaussian process regres- 
sion, kernel regression and even many heuristic homebrew 
algorithms. 

These algorithms all work on the idea of using the empir- 
ical training set to build up a full relationship between mag- 
nitudes and/or colours and the redshift. As each individual 
parameter (say the B — V colour) will have some spread with 
redshift, these give distributions or probabilistic values for 
the redshift, narrowed with each additional parameter. This 
process, in terms of artificial neural networks, is nicely de- 



scribed by Collister and Lahav ( 2004 1, who use this method 
in their publicly available photo-z code ANNz (described in 
the same paper). They also discuss the limitations and un- 
certainties that arise from this methodology. 

Machine learning algorithms (of which neural networks 
is one) are one of the strengths of the empirical method. 
These methods are able to determine the magnitude/colour 
and redshift correlations to a surprising degree, can handle 
the increasingly large training sets (i.e. SDSS) and return 
strong probabilistic estimates (i.e. well constrained uncer- 



tainties, see figure 16 1 on the redshifts (see|Ball et al. 2008a 



for a description of machine learning algorithms available 
and the strong photo-z constraints possible). In addition, 
machine learning algorithms are also able to handle the 
terascale datasets now available for photo-z determination 
rapidly, limited only by processor speed and algorithm effi- 
ciency ( |Ball et al.|2008b) . 

The additional benefit of the empirical method with ma- 
chine learning, now increasingly being used, is that the con- 
straining inputs for the photo-zs are not limited to the galax- 
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Fig. 15 Marginalized age distribution in four Virgo cluster el- 
liptical galaxies. Three models for the star formation history are 
considered as labeled (see text for details). The black lines corre- 
spond to the age extracted from [MgFe] plus either HyS (solid), Hy 
(dashed) or Hi5 (dotted), respectively. The gray line is the result 
from the fit to the SED. [Courtesy I. Ferreras] 




Fig. 16 Improvement in quasar redshifts enabled by the Bal l et al] 
12008b} data mining techniques, shown as spectroscopic versus 
photometric redshift for the SDSS. Left: A reproduction, using 
their framework, a typical result prior to their work (e.g. Weinstcin 
|eTaT{ 2004) Right: The result of using machine learning to assign 
probability density functions then taking the subset with a single 
peak in probability. Contours indicate the areal density of indi- 
vidual quasars (points) on the plot. Data from |B all et al. (2008b). 
[Courtesy N.Bali] 



ies SED. Suggested first by Koo ( 1999 1, properties such as 
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the the bulge-to-total flux ratio (e.g. 



surface brightness (e.g. Kurtz et al. 



Sarajedini et al.|1999) , 



20071, petrosian radii 



(e.g. Firth et al.|2003[), and concentration index(e.g. Collis- 



ter and Lahav 2004) have all been used in association with 



the magnitudes and colours to constrain the redshift, some 



codes even bringing many of these together (e.g. Wray and 
|Gunn|2008D . 

5.1.2 Template Fitting 

Unlike the empirical method, the template-based method is 
actually a form of SED fitting in the sense of this review (see 
e.g. |Koo||1985[ |Lanzetta et al.||1996| |Gwyn and Hartwick 
[T996}|Pello et al.|1996[|Sawicki et al.|T997] >. Put simply, this 
method involves building a library of observed (Coleman 



et al.|1 980 is a commonly used set) and/or model templates 



(such as Bruzual and Charlot|2003| l for many redshifts, and 
matching these to the observed SEDs to estimate the red- 
shift. As the templates are "full" SEDs or spectra, extrapola- 
tion with the template fitting method is trivial, if potentially 
incorrect. Thus template models are preferred when explor- 
ing new regimes in a survey, or with new surveys without a 
complementary large spectroscopic calibration set. A major 
additional benefit of the template method, especially with 
the theoretical templates, is that obtaining additional infor- 
mation, besides the redshift, on the physical properties of 
the objects is a built-in part of the process (as discussed in 
section |4~T|). Note though, that even purely empirical meth- 



ods can predict some of these properties if a suitable training 
set is available (see e.g. |Ball et al.|2 004). 

However, like empirical methods, template fitting suf- 
fers from several problems, the most important being mis- 
matches between the templates and the galaxies of the sur- 
vey. As discussed in section|2] model templates, while good, 
are not 100% accurate, and these template-galaxy colour 
mismatches can cause systematic errors in the redshift es- 
timation. The model SEDs are also affected by modifiers 
that are not directly connected with the templates such as 
the contribution of emission lines, reddening due to dust, 
and also AGN, which require very different templates (see 
e.g. |Pollettaetal.|2007| >. 



It is also important to make sure that the template set is 
complete, i.e. that the templates used represent all, or at least 
the majority, of the galaxies found in the survey (compare 



also Section 4.5.2 1. This is especially true when using em- 
pirical templates, as these are generally limited in number. 
Empirical templates are also often derived from local objects 
and may thus be intrinsically different from distant galaxies, 
which may be at different evolutionary stage. A large tem- 
plate set is also important to gauge problems with degen- 
eracies, i.e. where the template library can give two differ- 
ent redshifts for the same input colours. Another potential 
disadvantage of template fitting methods comes from their 



sensitivity to many other measurements to about the percent 
level, e.g., the bandpass profiles and photometric calibra- 
tions of the survey. 

For implementations of the template fitting, the method 
of maximum likelihood is predominant. This usually in- 
volves the comparison of the observed magnitudes with the 
magnitudes derived from the templates at different redshifts, 
finding the match that minimizes the x 2 (compare section 
4.5 1. What is returned is the best-fitting (minimum x 2 ) red- 
shift and template (or template +modifiers like dust attenu- 
ation). By itself this method does not give uncertainties in 
redshift, returning only the best fit. For estimations of the 
uncertainties in redshift, a typical process is to propagate 
through the photometric uncertainties, to determine what is 
the range of redshifts possible within these uncertainties. A 
good description of the template-fitting, maximum likeli- 
hood method can be found in the description of the publicly 



available photo-z code, hyperz in Bolzonella et al. (2000). 

As mentioned above, one of the issues of the templates 
is the possibility of template incompleteness, i.e. not hav- 
ing enough templates to describe the galaxies in the sam- 
ple. Having too many galaxies in the template library on 
the other hand can lead to colour-redshift degeneracies. One 
way to overcome these issues is through Bayesian inference: 



the inclusion of our prior knowledge (see Section 4.5 i, such 
as the evolution of the upper age limit with redshift, or the 
expected shape of the redshift distributions, or the expected 
evolution of the galaxy type fractions. As described in Sec- 
tion 4.5 this has the added benefit of returning a probability 



distribution function, and hence an estimate of the uncertain- 
ties and degeneracies. In some respects, by expecting the 
template library to fit all observed galaxies in a survey, the 
template method itself is already Bayesian. Such methods 
are used in the BPZ code of Bemtez (2000), who describes 
in this work the methodology of Bayesian inference with re- 
spect to photo-z, the use of priors and how this method is 
able to estimate the uncertainty of the resulting redshift. 

It should be noted that, while public, prepackaged 
codes might provide reasonable estimates for certain types 
of sources, no analyses should proceed without cross- 
validation and diagnostic plots. There are common prob- 
lems that appear in data sets and issues that need to be under- 



stood first, and worked around, if possible (see e.g. Mandel- 
baum et al.|2005 for a comparison of some public photo-z 
codes). Some further public photo-z codes include kphotoz 
( |Blanton et al.|2003| ), zebra ( |Feldmann et al.|2006| ) and he 
Phare ( |Arnouts et al.|1999||Tlbert et al.|2006||2009al ). 

5.2 Calibration and error budgets 

Redshift errors are ultimately data-driven: they typically 
scale with 1 +z given constant wavelength resolution of most 
filter sets; they also scale with photometric error in a tran- 
sition regime between ~ 2% and ~ 20%. Smaller errors 
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are often not exploited due to mismatches between data and 
model arising from data calibration and choice of templates, 
while large errors translate non-linearly into weak redshift 
constraints. If medium-band resolution is available, QSOs 
show strong emission lines and lead to deeper photo-z com- 
pleteness for QSOs than for galaxies. 

Photometric redshifts have limitations they share with 
spectroscopic ones, and some that are unique to them: as 
in spectroscopy, catastrophic outliers can result from the 
confusion of features, and completeness depends on SED 
type and magnitude. Two characteristic photo-z problems 
are mean biases in the redshift estimation and large and/or 
badly determined scatter in the redshift errors. Catastrophic 
outliers result from ambiguities in colour space: these are ei- 
ther apparent in the model and allow flagging objects as un- 
certain, or are not visible in the model but present in reality, 
in which case the large error is inevitable even for unflagged 
sources. Empirical models may be too small to show local 
ambiguities with large density ratios, and template models 
may lack some SEDs present in the real Universe. 

Remedies to these issues include adding more discrimi- 
nating data, improving the match between data and models 
as well as the model priors, and taking care with measur- 
ing the photometry and its errors correctly in the first place. 
Photo-z errors in broad-band surveys appear limited to a red- 
shift resolution near 0.02 X (1 + z), a result of limited spectral 
resolution and intrinsic variety in spectral properties. Trac- 
ing features with higher resolution increases redshift accu- 
racy all the way to actual spectroscopy. Future work among 
photo-z developers will likely focus on two areas: (i) Un- 
derstanding the diversity of codes and refining their per- 
formance; and (ii) Describing photo-z issues quantitatively 
such that requirements on performance and scientific value 
can be translated into requirements for photometric data, for 
the properties of the models and for the output of the codes. 

5.2.7 Template accuracy 

In general, template-based photo-z estimates depend sensi- 
tively on the set of templates in use. In particular, it has been 
found that better photo-z estimates can be achieved with an 
empirical set of templates (e.g. Coleman et al. 1980; Kin- 
ney et 5JT 996 1 rather than using stellar population synthe- 
sis models (SPSs; e.g. |Bruzual and Charlot|[2003| |Maras- 
ton||2005 see Section [2) directly. Yet the models are what 
are commonly used to compute stellar masses of galaxies. 
Since the use of these templates do not result in very good 
photometric redshifts, what is usually done, is to first derive 
photometric redshifts through empirical templates, and then 
estimate the stellar masses with the SSP templates. Obvi- 
ously this is not self-consistent. 

To investigate what causes the poorer photo-z estimates 
of synthetic templates, Oesch et al. (in prep.) used the pho- 
tometric data in 1 1 bands of the COSMOS survey ( |Scov- 
ille et al.||2007| >, together with redshifts of the zCOSMOS 



follow-up ( [Lilly et al.|2007) > and fit the data with SSP tem- 
plates. In the resulting rest-frame residuals they identified a 
remarkable feature around 3500 A, where the templates are 
too faint with respect to the photometric data, which can be 



seen in figure 17 The feature does not seem to be caused 
by nebular continuum or line emission, which they subse- 
quently added to the original SSP templates. Additionally, 
all types of galaxies suffer from the same problem, indepen- 
dent of their star-formation rate, mass, age, or dust content. 

Similar discrepancies have been found previously by 
[Wild et al.| ( p007] i; |Walcher et al.| ( |2008] i, who found a ~ 0.1 
mag offset in the D n (4000) index. As this spectral break 
is one of the main features in the spectrum of any galaxy, 
it is likely that the poor photo-z performance of synthetic 
templates is caused by this discrepancy. The cause of the 
discrepancy has been identified as a lack of coverage in the 
synthetic stellar libraries used for the models. It will thus 
be remedied in the next version of GALAXEV (G. Bruzual, 
priv. comm.). 

5.2.2 Spectroscopic Calibration of Photo-zs 

One of the strong benefits of the template method is that any 
spectroscopic subsample of a survey can be used to check 
the template-determined photo-zs. This can also be done for 
the empirical methods, yet for this a very large spectroscopic 
sample is necessary such that it can be divided into a large 
enough training set and testing sets. 

With the existence of a test spectroscopic sample, it is 
then possible to calibrate the template library, leading to a 
combined empirical-template method. This means to correct 
for errors in the photometric calibration or even the correc- 
tion of the templates themselves for example to allow for the 
evolution of galaxies with a small library, or to account for 
inaccurate models (see section [2]). Such calibration is typ- 
ically an iterative process, in which the photometry and/or 
template SEDs are modified to minimize the dispersion in 
the resulting photometric redshifts. 

The simplest kind of calibration involves adding small 
zero-point offsets to the photometry uniformly across the 
sample. This does not imply that the photometry is incor- 
rectly calibrated (though in practice the absolute calibration 
may well have small errors in the zero-point), but rather that 
there is often a mismatch between the real SEDs of galax- 
ies and the templates used to fit them. The calibration is 
meant to minimize those differences. Plotting color-color or 



color-redshift diagrams (figure 18 1 with the template SEDs 
overlaid will often indicate bulk offsets between the two. 

A more instructive approach, however, is to compute the 
residuals between the predicted magnitude of the best-fit 
template at the spectroscopic redshift and the observed mag- 
nitude (for more details, see |Brodwin et al.|2006a|b| l. These 
residuals can be plotted versus color or redshift for added 
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diagnostic power. In the example in Figure 19 there ap- 
pears to be an effective magnitude offset of * 0.3 mags in 
the //-band. 

Applying such effective zero-point adjustments in all 
bands in an iterative process minimizes the mismatch be- 
tween the data and the templates, and hence minimizes 
the resultant photometric redshift dispersion, as shown in 



Fig. 20 



Such calibration phases are used in the works of Brodwin 



et al. ( |2006a) > and as "template-optimization" in the codes 
zebra ( |Feldmann et al.|[2006) and Le Phare filbert et al.| 
2006 2009a) which use template fitting with Bayesian in- 
ferences and this calibration phase to give the most accurate 
photometric redshifts possible with the template approach. 

With the most accurate photometric redshifts possible, 
the template-fitting can then be used to estimate physical 
properties such as stellar masses, star-formation rates, etc. 
(see section[6]i. 

5.2.3 Signal-to-noise Effects 



Margoniner and Wittman ( 2008 1 have specifically investi- 
gated the impact of photometric signal-to-noise (SN) on the 
precision of photometric redshifts in multi-band imaging 
surveys. Using simulations of galaxy surveys with redshift 
distributions (peaking at z ~ 0.6) that mimics what is ex- 
pected for a deep (10-sigma R band = 24.5 magnitudes) 
imaging survey such as the Deep Lens Survey (Wittman 
et al.||2002 1 they investigate the effect of degrading the SN 
on the photometric redshifts determined by several publicly 
available codes (ANNz, BPZ, hyperz) 



Figure 21 shows the results of one set of their simula- 
tions for which they degraded the initially perfect photom- 
etry to successively lower SN. In these unrealistic simula- 
tions all galaxies have the same SN in all bands. The fig- 
ure shows the cumulative fraction of objects with 6 Z smaller 
than a given value as a function of S z . The left panel shows 
the cumulative fraction for all objects, while the right panel 
shows galaxies for which the BPZ photo-z quality param- 
eter, ODDS > 0.9. The number of galaxies in the right 
panel becomes successively smaller than the number in the 
left as the signal-to-noise decreases (64% of SN=250, and 
only 6.4% of SN=10 objects have ODDS > 0.9), but the 
accuracy of photo-zs is clearly better. 

The results of this work show (1) the need to include 
realistic photometric errors when forecasting photo-zs per- 
formance; (2) that estimating photo-zs performance from 
higher SN spectroscopic objects will lead to overly opti- 
mistic results. 

5.3 A unified framework 

The field of photometric redshifts and the estimation of other 
physical properties has been very pragmatic. Its develop- 
ment, since the first attempts, has been incremental in the 



sense that most studies focused on refining components but 
staying within the concepts of the original ideas of the two 
classes. Empirical and template-fitting approaches today 
follow very separate routes, with these classes of methods 
even use different sets of measurements. Only the semi- 
empirical approach of zero-point calibration comes close to 
linking the two approaches. However, a recent study by 



Budavari ( 2009a i, tries to understand this separation and 
possibly bring these methods together by devising a unified 
framework for a rigorous solution based on first principles 
and Bayesian statistics. 

This work starts with a minimal set of requirements: a 
training set with some photometric observables x and spec- 
troscopic measurements and a query or test set with some 
potentially different set of observables y. The link between 
these is a model M that provides the mapping between x and 
y, the probability density p(x\y, M). This is more than just 
the usual conversion formula between photometric systems 
because it also incorporates the uncertainties. 

The empirical relation of x-% is often assumed to be a 
function. A better approach is to leave it general by mea- 
suring the conditional density function. The simplest way 
is to estimate the relation by the densities on the train- 
ing set as p(t;\x) = p(£;,x)/p(x). The final result is just 
a convolution of the mapping and the measured relation: 
p(g\y,M) = Jdx p(g\x)p(x\y,M). In figure 22 we show the 
results from Budavari (2009a), where he plots the empirical 
relation (dotted blue line) and the final probability density 
(solid red) for a handful of SDSS galaxies. The top panels 
show intrinsically red galaxies, whose constraints are rea- 
sonably tight out to the highest redshifts. Blue galaxies in 
the bottom panels however get worse with the distance as 
expected. 

The aforementioned application follows a minimalist em- 
pirical approach but already goes beyond traditional meth- 
ods. Template-fitting is in the other extreme of the frame- 
work where the training set is generated from the model us- 
ing some grid. Without errors on the templates, the equa- 
tions reduce to the usual maximum likelihood estimation 
that is currently used by most codes. A straight forward 
extension Budavari (2009a) suggests is to include more re- 



alistic errors for the templates. Similarly one can develop 
more sophisticated predictors that leverage existing training 
sets and spectrum models at the same time. 

5.4 The State of Photometric Redshifts 

Generally there is an overall agreement in most aspects of 
photometric redshift methodologies, and even technicalities. 
However there is a need for standardized quality measures 
and testing procedures. It is important to analyze to perfor- 
mance of each model spectrum as a function of the redshift. 
This is best done by plotting the difference of the observed 
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magnitudes and template-based ones. These figure can pin- 
point problems with the spectra and even zero-points. In- 
trinsically these are the quantities used inside the template- 
optimization procedures, e.g., in |Budavari et al. (2000! and 
|Feldmann etaT1 ( |2006l l. 

In SED fitting, interpolation between templates is often 
used, which can be done linearly or logarithmically. The lat- 
ter has the advantage of being independent of the normaliza- 
tion of the spectra. Yet, most codes appear to use linear in- 
terpolation without a careful normalization. This might ex- 
plain some of the discrepancies among similar codes found 
by the Photo-Z Accuracy Testing (PHAT) project. |^| 

The determination of the quality of the estimates is also 
a crucial topic. There is need for different measures that 
can describe the scatter of the points without being domi- 
nated by outliers and that can estimate the fraction of catas- 
trophic failures. It is also recommended to characterize the 
accuracy of the estimates by a robust M-scale instead of the 
RMS; a measure that is simple to calculate, yet, not sensitive 
to outliers. Another aspect of this is the study of selection 
criteria that is often neglected. Certain projects are not con- 
cerned with incomplete samples as long as the precision of 



the ones provided is good (e.g., weak lensing, Mandelbaum 
et al.|2005) , while others, such as galaxy clustering, might 



require an unbiased selection. Therefore, it is perceived that 
studies using methods with any quality flags or quantities 
should provide details of their selection effects. 

A common theme for future goals in most photomet- 
ric redshift works appears to be more detailed probabilis- 
tic analyses, with the need for probability density functions. 
Priors used in most bayesian analyses seem to be generally 
accepted in the photo-z community. With such consensus 
amongst photometric redshifts obtained, the focus of work 
now is shifting from the estimation of "just" the redshifts 
to simultaneously constraining physical parameters and the 
redshift in a consistent way. 



6 Results of SED Fitting: Physical Properties 

SED fitting is a very versatile tool. From a rough estima- 
tion of the stellar masses of distant galaxies to the search 
for small subpopulations of stars in high S/N spectra, it can 
be applied to a large part of the problems in galaxy evolu- 
tion. This is the strength and the weakness of the SED fitting 
technique: it does it all at once. 

We highlight here a few significant results. The intention 
is not to be complete or to mention the work that has been 
most in view, but rather to highlight the diversity of ques- 
tions that can be adressed from fitting the integrated SEDs of 



stellar populations. Particular importance has been given to 
supply cautionary remarks, as it is easy to overinterpret the 
significance of the derived properties, in view of the com- 
plexity of the physical mechanisms and our frequent lack of 
detailed understanding. 

6.1 Stars 

6.1.1 Stellar masses 

Stellar masses are computed by multiplying a mass-to-light 
ratio M/L with a luminosity L. While the uncertainties on L 
depend on the quality of the data, the estimate over M/L and 
its associated uncertainties depend mostly on the care taken 
with SED fitting. It is a good idea to search for a reference 
band that minimizes the effects of M/L variations due to 
stellar populations (age, metallicity, chemical abundances) 
and due to dust absorption. While the common notion that 
the NIR (e.g. the H band at 1.65 pun) is close to ideal is 
correct in some cases, problems can arise because of ther- 
mally pulsing asymptotic giant branch stars (discussed in 
section 2. 1 .3 1 if young ages (<2 Gyr) are well represented in 
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the SED. That reliable M/L from SED fitting cannot be dis- 
pensed with is evident when looking at IRAC 3.6jum data of 
nearby galaxies, where star formation regions are evidently 
prominent. 

Stellar Mass Maps of Resolved Galaxies In the work of 
Zibetti et al. (2009) a "data-cube" approach is introduced 
to investigate the SEDs of nearby, resolved galaxies, aimed 
at preserving the maximum spatially-resolved information. 
One feature of the approach is that it allows to compare 
"global" quantities, which are notoriously difficult to deter- 
mine, with the integral over the local quantities, a useful test 
of how meaningful global quantities can be. A large part 
of the effort concentrated on developing a reliable method 
to obtain stellar surface mass density maps from a mini- 
mum set of broad-band observations. This method relies on 
Bayesian inference, as discussed in Section [45] 

The effective mass-to-light ratio correlates with optical / 
near-infrared colors (e.g., Bell and de Jong 2001), so M/L 
can be expressed as a function of color(s). A better estimate 
is obtained if M/L is mapped as a function of two colors, 
instead of one. The colors adopted are g - i and i - H. Their 
large wavelength separation allows to robustly describe the 
shape of the SED over the entire optical-near-infrared range, 
in a way that as insensitive as possible to photometric and 
modeling uncertainties. 

To study the dependence of M/L on (g—i,i—H) the au- 
thors use a Monte Carlo library of 50,000 models cre- 
ated from the 2007 version of BC03, which include also 
dust in different amounts and spatial distributions.. The 
(g—i,i—H) space is binned in cells of 0.05 mag x 0.05 mag 
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and marginalized over M/L in each cell. The median is cho- 
sen as the fiducial M/L at each position of the color-color 
space. A look-up table is created to derive M/Lh as a func- 
tion of (g—i,i—H). Figure 23 illustrates that the information 
in the second color improves the M/L determinations sys- 
tematically by +0.3-0.4 dex. 

The method described above is applied to each pixel of 
the image of a galaxy, where "pixel" implies the pixel that 
results after matching the images in the three bands to a 
common resolution. In order to provide sensible results, it is 
crucial that color measurements do not exceed 0. 1 mag er- 
rors, which requires S/N-20. The results for M 100 (NGC 



4321) are displayed in Figure 24 



An important result appears from the comparison be- 
tween total stellar mass as obtained by integrating the stel- 



lar mass surface density maps (Figure 24 1) and the one ob- 
tained using global colors to estimate the "average" M/L 
ratio to be multiplied by L(H). This second method, the one 
commonly adopted in extragalactic studies, agrees within 
~10% of the mass map integral only when the color dis- 
tribution is quite homogeneous, i.e., for early type galax- 
ies. When substantial color inhomogeneities and especially 
heavily obscured regions are present within a galaxy, using 
global colors and fluxes can lead to underestimates of the 
actual stellar mass of a galaxy by up to 60%. This can be 
understood if dust-obscured regions can contribute a signif- 
icant amount of mass but are heavily under-represented in 
the global color, which is flux weighted and hence biased to- 
ward the brightest (and bluest) parts of the galaxy. While the 
pixel-by-pixel M/L gives the correct mass weight to these 
obscured regions, the globally computed M/L severely un- 
derestimates their mass contribution. 

Stellar mass functions One of the holy grails of current ob- 
servational efforts in galaxy evolution studies is a consistent 
picture of the build-up of stellar mass over the age of the 
universe. An important constraint on this is the stellar mass 
function, or its integral, the stellar mass density. A compre- 
hensive discussion of these results would warrant a review 
of its own. Suffice it here to point out that not only the lo- 
cal mass function has been measured with great precision 
(e.g. Bell et al.|2003a i, but these results have also been ex- 
tended to redshifts of 1 (Pozzetti et al. 2007 Bundy et al. 
2006). At redshifts above 1.2 an observed-frame optical se- 
lection corresponds to a rest-frame UV selection, subject to 
large biases. These therefore have to be circumvented by a 
K-band selection (e.g. Cirasuolo et al.|2 007) or, better, by a 
selection at 3.6/im ( |Arnouts et al.||2007| [Tlbert et al.||2010| 
e.g.). For observational and conceptual reasons (detailed 



in Sections 4.5.2 and 4.6 1, determining stellar masses and 



therefore mass function at redshifts higher than z=2 is very 
difficult. Most authors thus prefer to restrict themselves to 
luminosity functions instead (see e.g. Cirasuolo et al.|20T0 



for just one very recent example), thus leaving the burden of 
transforming luminosities to stellar masses to the interpret- 
ing model (but see e.g. Kajisawa e t al.|2009 i. 



Stellar masses of high redshift galaxies Many authors have 
used some kind of Bayesian inference-based method (Sec- 



tion 4.5 1 to determine stellar masses for high redshift galax- 



ies (e.g. |Sawicki and Yee 1998 Pa povich et al.|2001j|Forster 
Schreiber et al.|2004) >. There is good hope that these mass 
estimates are reasonably good ( |Drory e t al. 2004), despite 
important caveats on the methodology that become more 
important with improving data quality (see e.g. Section 
4.1). A recent result has been the discovery and study of 
high redshift galaxies with high stellar masses and low star 
formation rates (early type galaxies, ETGs). In the follow- 
ing we describe only one "family" of papers, as presented 
at the workshop, but see Cimatti (2009 ) for a review. Mas- 
sive ETGs are the first objects to populate the red sequence 
(see e.g. |Kriek et al. 2008b). Objects in the redshift range 
2-3 can be identified by multi-band photometry (e.g. 



van 



Dokk um et al.|2006 >. For determining the physical proper- 
ties however, significant uncertainties are due to photomet- 
ric redshift determinations. For example [Kriek et aT| ([2008 ) 
conclude that while stellar masses are reasonably robust to 
small errors arising from photometric redshifts, the actual 
star formation history is generally very poorly constrained 
with broad band data alone. 

The obvious next step is thus to analyze these galaxies 
with spectroscopy, despite this being an expensive under- 



taking in terms of telescope time. When doing this, Kriek 



et al. ( 2008 1 also go further in blurring the limits between 
spectroscopy and photometry by binning their "low" S/N 
spectra into bins of 400 A. While the information content 
remains unchanged, this certainly leads to improvements in 
presentability and fitting speed. Kriek et al. (2009a) then 



show explicitly that provided enough photons can be assem- 
bled through either exposure time or telescope size, the spec- 
tra of galaxies at redshifts 2-3 are amenable to the same kind 
of analysis as in the local universe. The upshot of these stud- 
ies is that massive, compact ETGs with very little residual 
star formation are in place already at redshifts between 2 and 
3. 

Despite these successes, the study of Muzzin et al. 
( 2009a b) confirms that even using spectroscopic data, 
model uncertainties mean that SED-derived stellar masses 
are affected with uncertainties of factors 2-3 at these red- 
shifts. For further discussion on stellar mass determinations 
the reader is also referred to the review by de Jong & Bell 
(in prep.). 

6.7.2 Deriving SFHsfrom spectroscopy 



Comparing observations to semi-analytic models Trager 



and Somerville ( 2009 ) extend the semi-analytic model of 
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Somerville et al. 



(2008) to predict the line absorption 



strengths (Section 4.2 1 of the resulting galaxies. This allows 
them to use the same analysis tools that would be used in the 
analysis of the measured line strengths of an observed sam- 
ple of galaxies on objects with known properties, in partic- 
ular star formation histories. They select in particular early 
type galaxies from the mock catalogues they produce and 
compute the index strengths of the resulting spectra. These 
index strengths can then be plotted on the same plots as real 
data. They come to the sobering conclusion that while the 
sample of Trager et al. ( 2008) is of sufficient quality do do a 
meaningful comparison, it remains too small. On the other 



single tracers (see Section 4.5.3 for a comparison). Walcher 
|et al.| ( |2008| l have used stellar masses and star formation 
rates consistently derived from the same SED fit to compare 
the predicted evolution of the stellar mass function to the ob- 
served one. The main result is that while stars form in blue 
cloud galaxies, most of the growth of the stellar mass func- 
tion occurs in quiescent galaxies, in agreement with studies 



hand large samples of galaxies, as the one of Moore et al 



based on different tracers of star formation (e.g. Bell et al. 
2007| l. From comparison with merger studies in the same 
field, jWalcher efaX] ( |2008| > conclude that about half of the 
mass growth on the red sequence comes from major mergers 
and half from minor mergers. 

Salim et al.|(j2007 ) have compared their SED-fitting SFRs 



( 2002 1, still lack the required precision. 



The archeology of the universe The database of the SDSS 
spectra has been used to derive the SFH of galaxies from 
their current spectra (e.g. Heavens et al. 2004 see also 



Section 4.4 1, a procedure sometimes called "unlocking the 
or simply "astro-archeology". A recent up- 
Tojeiro et al. ( 2007), who 



fossil record 
date on this has been presented in 



applied VESPA (see Section 4.4 



to the SDSS sample of 



to SFRs determined from emission lines. They find that 
some galaxies with no detected emission lines nevertheless 
have substantial SED-based SFRs. They attribute this result 
to "recent" star formation, i.e. stars that formed long enough 
ago that their emission lines already vanished, but still re- 
cently enough to be revealed in the galaxy SED. Recent HST 
imaging in the UV which clearly shows SF structures seems 
to confirm this (Salim & Rich 2010, ApJ submitted). 



spectra and derived a catalogue that was made available to 6 ] j identifying and studying outliers 



the community at http://www-wfau.roe.ac.uk/vespa/ It pro 



vides detailed star formation histories and other parameters 
for SDSS's latest data release (DR7) of the Main Galaxy 
Sample and the Luminous Red Galaxies sample. Details 
of the catalogue, including description, basic properties and 
example queries can be found in ( |Tojeiro et al.|2 009). 



Combining spectroscopic and dynamic ages The use of 
spectroscopy in combination with dynamical arguments to 
understand the evolution of a single object was presented in 
Pappalardo eTaT1 ( [2010) l. The galaxy NGC4388 is a member 
of the Virgo cluster and sports a huge trail of HI gas (Oost-| 



erloo and van Gorkom 2005). It represent an ideal study 
case for the effects of ram stripping on a galaxy moving in 
the intra-cluster medium (Vollme r and Huc htmeier 2003}. 
The stripping age has been estimated to be of order 200 
Myr from dynamical arguments. Using VLT/ FORS spec- 
troscopy of the outer and inner regions of the galaxy and 
the STECKMAP program ( |Ocvirk et al.|2006) , |Pappalardo 
et al. ( 2010) 1 were able to show that, while the inner region of 
the galaxy is of solar metallicity and has continued forming 
stars to the present day, the outer regions of the galaxy have 
sub-solar metallicity and have stopped forming stars roughly 
200 Myrs ago, in accordance with the dynamical estimate. 
Single cases like this can thus help identify the processes 
and timescale associated with shutting down the star forma- 
tion in galaxies, one of the most profound changes a galaxy 
can experience. 

Star formation rates Star formation rates from SED fitting 
have been little used, with most authors preferring to rely on 



This is an underexplored use of SED fitting, in the opinion 
of the authors. One example, objects with differing SFR 
measurements from emission lines and SED fitting has been 
covered in the last section. 

Finding Wolf-Rayet stars The availability of large 
databases of spectra, such as from the SDSS, and of 
accurate stellar population model predictions enables the 
search for rare objects or systematic deviations that are 
not predicted by the model itself. An example of this are 
Wolf-Rayet stars, evolved, massive stars with characteristic 
features. These have ages between 2 x 10 6 and 5 x 10 6 
years, and are thus a transient feature of galaxy spectra. 
They are useful as a tracer of recent star formation history 
as well as the metallicity of their host galaxies. As shown in 
Brinchmann et al. ( 2008| l, a systematic search in the SDSS 
database yields a sample several times larger than previous 
serendipitous searches. The essential ingredient of such a 
search is the accuracy of the stellar population model that 



allows an inversion technique (Section 4.4 1 to be applied 
on a large sample. Indeed, either a smooth correction or 
residual features from inaccurate models would severely 
impair the identification of the specific features. As an 
example Brinchmann et al.| (|2008|) show that the use of 



the Bruzual and Chariot ( 2003) models produced a large 
number of false positives, while an updated version of the 
same model using different stellar spectra (CB09) provides 
much better fits. 
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6.2 Dust 

Dust cannot be ignored when fitting a galaxy's SED, as 
shown by the cosmic infrared background, which has com- 
parable power to the distinct peak of the cosmic UV-optical 
background (Hauser and Dwek 2001). The relative strength 
of the cosmic background in the infrared suggests a signif- 
icant processing of the galactic stellar light over the age of 
the universe. This processing must have also been more sig- 
nificant with increasing redshift as the percentage of stellar 
light re-radiated by dust is only ~ 30% locally (Popescu and 



The 217 5 A feature The 2 175 A feature has been associated 



with small carbonaceous grains in the ISM (Mathis et al. 



|Tuffs|2002l >, as supported by the increasing number density 
of luminous IR galaxies up to z ~ 1.3 (Magnelli et al. 2009 ). 

As discussed in section[2] the absorption and emission of 
light by dust are generally treated as separate processes in 
modelling, and this is similarly true in SED fitting. 

6.2.1 Attenuation by Dust 



Dust between the observer and the individual stars of a 
galaxy acts to extinguish and redden the light from those 
stars. When the stars in our own Galaxy were examined it 
was found that a simple relation with wavelength was able 
to describe the extinction and reddening by dust for a wide 
range of galactic envir onments, with the onl y strong feature 
occurring at -2175A (Cardelli et al.||l989i. A similar but 



steeper extinction law was found for the Magell anic Clouds, 
with weaker or non-existent feature at 2 175 A (Gordon and 



Clayton|19 98 1 Misse lt et al.|1999| >. It is these extinction laws 
that have given rise to the contemporary model of dust in the 
ISM (i.e. Math is et al.|1977 1, and the understanding that the 
dust composition between the Milky Way and Magellanic 
clouds is different. 

Yet when integrated over the whole of a galaxy the situ- 
ation becomes complex, with the geometry of the stars and 
dust strongly affecting the resulting spectrum. The effects of 
varying amounts of extinction of the different stellar popula- 
tions due to the spatial distribution of stars and clumpy dust, 
and the scattering of blue stellar light into our line of sight 
act to flatten the effects of dust on the spectrum, creating an 
attenuation law, where the amount of reddening with extinc- 



tion is less (or 'greyer') than we observe locally (Witt et al. 



1992). This was exactly what was found in starburst galax- 



ies by |Calzetti et al] ( [1994) ), and |Charlot and Fall] pOOO i 
found that a simple screen effective attenuation (i.e. a screen 
of dust between the galaxy and observer) with a power-law 
relation, Tism k /l~ ' 7 , was able to account reasonably well 
for the diffuse ISM attenuation in galaxy observations. It is 
this complexity that makes disentangling the effects of ge- 
ometry and differing dust difficult, and thus the extraction 
of physical dust properties from galaxy SEDs problematic. 
There are two areas where some progress has been made. 



1977| ), and is observed in both the Milky way, M33 (Gor- 



don et al. 1999), and (weakly) in the LMC, but is almost 
non-existent in the SMC. This feature is not observed in the 
attenuation law of starburst galaxies ( |Calzetti et al.|[T994| . 
Whether this lack is due to the clumpy geometry of dust and 
stars (Fisc hera et al.|2003] ) or is actually indicative of SMC- 
like dust in starburst galaxies (Gordon et al. 1999[ ) is still 
under debate, yet this feature is generally not needed to fit 
the attenuation of galaxies. In QSOs, which, being domi- 
nated by a nuclear source, are closer to the galactic extinc- 
tion situation, an average attenuation curve does not show 
this feature, suggesting processing of the ISM in these active 
objects ( Czerny et al.||2004] >. However, in a few non-local 
galaxies where direct extinction lines of sight are available, 
this feature has been observed, suggesting it may be more 
common than the attenuation curves of local galaxies sug- 
gest ( |Wang et al.][2004t |Eliasd6ttir eraLl[2009) ). At higher 
redshifts, where UV spectra are more commonly observed, 
recent studies find evidence for the existence of the 2 175 A 
bump (e.g. |Noll et al.|20"07t|Noterdaeme et al.|2009) . 



Young versus old attenuation One important progress 
made in the treatment of galaxy attenuation is the realisa- 
tion that the effective attenuation of a galaxy is dependent 



upon its star formation history. Calzetti ( 1997 1 found that in 
starburst galaxies the effective attenuation of the stellar con- 
tinuum was less than that suffered by the nebula emission, 
as measured through emission lines. This clear indication of 
the dumpiness of the dust in galaxies has been interpreted 
as an indication of differential attenuation of different stel- 
lar populations, with young stars, and their associated ion- 
ized nebula, strongly attenuated by the clouds from which 
the stars formed, while older stars have evolved out of their 
'birth clouds' either through cloud or stellar dispersion, and 
are only attenuated by the diffuse ISM dust, which acts on 



both the young and old stars (see e.g. Chariot and Fall 2000). 
Exactly what is the clearing time of these clouds and the 
differential attenuation is still uncertain, and may be galaxy 
specific, but this forms a basis for current galaxy SED mod- 
els as discussed in section [2] 

6.2.2 Dust Emission 



Extracting physical properties from dust emission in the IR 
is difficult for both theoretical and observational reasons: 
excluding the mid-IR, there are no observed dust features, 
most being washed out due to the broad shape of the black- 
body emission; the IR suffers from strong observational con- 
straints, with most data coming from space- and balloon- 
based observations; associated with this is the, until recently, 
limited sensitivity and spatial resolution and in the far-IR, at 
wavelengths > 100/rni, the sparsity of data. 



38 



With ISO and, especially, Spitzer space telescopes this 
situation has recently improved, and will improve more so 
in the near future with the recent launch of Herschel and 
ALMA beginning to take form. So we touch upon here some 
of the galaxy physical properties that have been determined 
from the dust IR emission. 

PAH emission in the Mid-IR As mentioned in section [2] 
the 5-20//m mid-infrared spectrum of galaxies is gener- 
ally dominated by broad emission features arising from 
large molecules, polycyclic aromatic hydrocarbons (see e.g. 
Smith et al. 2007 ). Underlying these features is the stellar 
continuum at short wavelengths and hot dust emission. Con- 
fusing the interpretation of the emission features are strong 
ionic emission lines arising from species such as Ne + and 
Ne ++ and strong, broad absorption features from silicate 
grains at 9.8yimi and 18/im. 

A recent tool, PAHFIT, has been developed to decom- 
pose the mid-infrared spectra into its stellar, PAH, dust 
continuum, and line emission constituents, using functional 
forms and templates for the features in this wavelength range 
( Smith et al. 2007| l. An example of this can be seen in figure 



25 The PAH feature luminosity has been used as star forma- 



tion rate tracers (see section 6.2.4 1, and the relative strength 
of these features to the continuum have been found to be 
strongly linked to the presence of AGN (see e.g.|Spoon et al. 



2007 [ and below), and to the gas phase metallicity (see e.g. 



Smith et al.|2007 1. The relative strengths of these features 



can also be used to diagnose the mean size and ionization 
state of the PAHs, which is related to the average radiation 
field and dust size distribution (D raine and Li|2007a| >. 

Diagnosing the energy source in ULIRGs Due to the high 
obscuration by dust in IR bright galaxies, especially ultra 
luminous IR galaxies (ULIRGs), diagnosing the dominant 
heating source is problematic. The diagram put forward 
by Spoon et al. (2007]) helps resolve this issue by using the 
strength of the strong silicate absorption feature that is de- 
termined from fitting the mid-IR SED (as discussed above 



in Section 6.2.2 1 in association with the equivalent width 
of the PAH features. This diagram cleanly separates dif- 
ferent classes of ULIRGs, from obvious Seyfert galaxies, 
strongly starbursting galaxies, and to deeply buried AGN 
ULIRGs and represents one of the strengths of IR SED fit- 
ting, extracting information from objects which are heavily 
obscured at shorter wavelengths. 

Dust masses One of the more important properties ob- 
tained by fitting the IR SED is the dust mass. Through fit- 
ting of the far-IR SED the temperature(s) and the relative 
contributions of the different dust components that make up 
the SED can be constrained. Then, using knowledge of the 
emissivity per unit mass of dust, the total dust mass (M<j) can 



be determined, using an equation such as (based on Dunne 
land Eales|20"0l] ); 



Mi = 



L&50 

*d(850) 



E 



B(850, T k ) 



(5) 



with Ls50 the 850jt/m luminosity, and Nk and Tk the relative 
contribution and temperature of dust component k. The sum 
of dust components is usually limited (< 3) by the sparse 
observational points at long wavelengths, but can also be 
represented by an integral of temperatures, parametrized by 
the strength of the heating radiation field (such as used by, 
e.g. Dale and Helou|200"2} |Draine and Li|2007a| >. Kd(850) is 
the dust mass opacity coefficient, taken to be 0.077 m 2 kg~' 
by Dunne et al. ( 2000| ); Dunne and Eales (2001), an interme- 
diate value between graphite and silicate. It is generally with 
this parameter that most of the uncertainties in determining 
dust masses remain. 

Longer wavelength fluxes (> 300pm, such as 850jt/m) 
are preferable to shorter wavelengths when determining 
dust masses as these sample the Rayleigh- Jeans part of the 
Planck curve, where the flux is least sensitive to temperature. 
Longer wavelengths are also more sensitive to the mass of 
the emitting material, as they are sensitive to cold dust as 
well as warm. 

Clear examples of fitting the far IR SED using simple, 
emissivity-modified black-bodies and determining the total 
dust masses can be found in |Dunne et al.| ( |2000) ) jDunne and| 
Eales ( 2001[ ) and more recently in Clements et al.| (2009 ) 



(see also da Cunha et al.||2010| Section 63 below). These 
works detail nicely the pertinent issues with both the data 
and fitting the far-IR SEDs. One of the best examples of de- 
termining the total dust mass, as well as other parameters, 



using the full IR SED was done by Draine et al. (2007b). 
Their physically based SED models (described in detail in 



Draine et al. 2007b ) were fitted to the far-IR SEDs of galax- 



ies from the SINGS sample, and gave determinations of the 
total dust mass, PAH fraction and information on the inter- 
stellar radiation field heating the dust. They found that dust 
in spiral galaxies resembled that found in the local Milky 
Way ISM, with similar dust-to-gas ratios, and that gener- 
ally it is the diffuse ISM that dominates the total IR power, 
excluding strong starbursting systems. These results thus 
confirmed the earlier ISO discoveries (see the review by 
Sauvage et al.|2 005 ). Note also that even earlier detailed ra- 



diative transfer modelling of individual galaxies had pointed 
to the dominance of the diffuse component ( |Popescu et al. 
|2000l >. 

Sub-mm excess emission SED fitting can not only return 
physical properties, but can also indicate where our current 
knowledge is failing. As mentioned above, the long wave- 
length dust emission is a good handle for the total dust mass. 



39 



However, when fitting the IR SED of several dwarf galax- 
ies it has been found that the sub-mm flux is in excess to a 
standard cool dust-body emission, requiring additional dust 
components at a unreasonably low temperatures (< 7 K) to 
fit the SED ( |Lisenfeld et aLpOO^ |Israel et aLpOTO] e.g.). 
While very cold large grains could be one possible cause, 
other suggestions have included small stochastic grains that 



spend most of their time at cold temperatures (Lisenfeld 
et al.|2 002 ), rotating dust grains (Israel et al. 2010| l, or some 



modification of the dust emissivity at these wavelengths or 



temperatures (Draine and Lee 1984 Weingartner and Draine 
|2001| l. Either way until this issue is resolved on the cause 
of this excess, the dust mass of these dwarf galaxies such 
as NGC 1569 will have large uncertainties. It is hoped that 
telescopes such as Herschel and ALMA may find more of 
these objects in the near future and help find the cause of 
this excess emission. 

6.2.3 Dust in the UV to IR 

The infrared-to-ultraviolet ratio is a coarse measure of dust 
extinction in the ultraviolet, and thus should be related to 
the amount of reddening in ultraviolet spectra. Indeed, star- 
bursting galaxies follow a tight correlation between the ratio 
of infrared-to-ultraviolet emission and the ultraviolet spec- 
tral slope (e.g. Calzetti 1997 Meurer et al. 1999| l. Compared 
to the relation defined by starbursts, normal star-forming 
galaxies are offset to redder ultraviolet spectral slopes, ex- 
hibit lower infrared-to-ultraviolet ratios, and show signif- 
icantly larger scatter ( |Kong et aL"1|2004[ |Buat et al.1|2005[ 
Burgarella eFaT1[2005] |Seibert et aE1[2005] |Cortese et a~ 



2006 , Boissier et al. 2007; G il de Paz et al.|2 007; Dale et al 



2007 ). Offsets from the locus formed by starbursting and 



normal star-forming galaxies can be particularly pronounced 
for systems lacking significant current star formation, such 
as elliptical galaxies, systems for which the luminosity is 
more dominated by a passively evolving older, redder stellar 
population. 

Using a sample of 1000 galaxies with spectroscopy from 
the SDSS and homogeneous photometric coverage from the 
UV to 24yum from SDSS and the Galex and Spitzer satel- 



lites, Johnson et al. ( 2007a i found that the sample galax- 
ies span a plane in the three-dimensional space of NUV- 



3.6/im colour, D„(4000) index (as defined by Balogh et al. 
|1999| l, and infrared excess, IRX (=Lm?/Lf j/y). The three- 
dimensional relation can be expressed in terms of empirical 
functions, where IRX is a function of NUV-3.6pm (or more 
weakly with other colours) and D„(4000). They suggest that 
this relation can be explained primarily through SFH and 
dust attenuation, with both acting to steepen the optical-UV 
slope (as measured by the NUV-3.6/im color), but only at- 



tenuation increasing the IR flux and hence IRX (Johnson 
|et al.|2006] >. 



A similar analysis was presented at the workshop by D. 

which con- 



Dale using the LVL survey (see Section 3.3.1 



sists of a statistically complete set of star-forming galaxies, 
nearly two-thirds of which are dwarf/irregular systems. Fig- 
ure [26] shows the ratio of far-ultraviolet-to-near-infrared lu- 
minosity as a function of the (perpendicular) distance from 



the starburst curve (e.g. |Calzettip 997 l Me urer et al. || 1 999 1 
for the LVL galaxies, with the far-ultraviolet emission is 
corrected for attenuation using the infrared-to-ultraviolet- 
based recipe formulated in |Buat et al. ( 2005| l. By correct- 
ing for dust, the FUV/3.6 fim ratio measures only the ra- 
tio of past-to-present star formation, sometimes referred to 



as the birthrate parameter (see also, for example, Boselli 
et al.||200"Tj Cortese et al.||2006| l. This ratio represents the 
birthrate parameter since the far-ultraviolet traces star for- 
mation over 100 Myr timescales whereas the near-infrared 
probes the total stellar mass built up over much longer 
timescales. This plot is as such a compression of the plane 



discussed by |Johnson et a l. (2007a), and shows a clear trend, 
with lower birthrate systems exhibiting larger distances from 
the starburst trend, consistent with the study of Kon g et aT| 
(2004). To further quantify this, theoretical models with so- 
lar metallicity, 1 M Q yr _1 continuous star formation curves 
assuming a double power law initial mass function, with 
Q'i.imf = 1-3 for 0.1 < m/M Q < 0.5 and q^.imf = 2.3 for 
0.5 < m/M G < 100 were run (Vazquez and Leitherer|2005 



and were matched with their determined FUV/3.6 /im ratio 
on the right axis, demonstrating that those with the oldest 
SFH (i.e. lowest birthrate parameters) lie furthest from the 
theoretical starburst curve. 

6.2.4 Star Formation Rate from the IR 

One of the most commonly extracted galaxy properties from 
the IR is the star formation rate. As discussed above, re- 
cently formed stellar populations tend to be more obscured 
than older stellar populations. They are also more luminous 
and emit more in the ultraviolet where dust opacity peaks, 
and thus dust emission is in principle a good tracer of recent 
star formation, assuming a simple calorimetric situation. It 



is these assumptions that lead to the widely-used Kennicutt 
( 1998 ) relation between total IR luminosity (8-1000//m) and 
star formation rate. 

Of course the situation is more complex than this, as dis- 
cussed in the same work. A range of ages contribute to dust 



heating ( |Kennicutt et al. 2009) 1, and star forming regions 
in galaxies suffer a range of obscurations, from totally ob- 
scured (ULIRGs) to unobscured (blue compact dwarfs). It 
is for this reason that this relation has been re-examined and 
empirically calibrated with new data from Spitzer. In gen- 
eral, all recent studies have found that the IR can be used 
as a SFR indicator, albeit without a direct one-to-one con- 
nection. Complications arise due to the correlation of SFR, 
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luminosity, and galaxy gas and dust masses, and possible 
non-linearities due to metallicity effects ( Wu et al.|2005 1. 
More specifically, |Calzetti et al. ( 2007 1 using spatially re- 



solved observations, found strong correlations between the 
24/mi and Pa luminosity densities (a proxy for SFR, as- 
suming little attenuation in the near-IR), and a correlation 
between the 8//m and Pa luminosity densities, though this 
failed at low metallicities. Using these, they were able to 
create new calibrations for SFR versus 24 fim luminosity, 
and SFR versus 24 fim and observed Ha luminosities, with 
the latter relation accounting better for the escaping radia- 
tion not accounted for by the dust emission. Rieke et al. 



(2009) took this further, showing that for higher IR lumi- 
nosity objects, Pa was no longer a good tracer for SFR as 
even it was obscured, and gave their own calibrations for 
SFR with the IR luminosities. 

On galactic scales, Zhu et al. ( 2008 ) showed with a larger 
galaxy sample from SWIRE that the Calzetti et ah] p.007\ 
relations between 24/im luminosity and extinction corrected 
Ha luminosity hold, and thus L(24/vm) is a good SFR indi- 
cator. In addition they also showed that 24yum is well cor- 
related with 70/mi and total IR luminosities, indicating that 
these too can be used as SFR indicators, albeit with larger 
scatter. 

The 8/im (and other PAH bands) and longer wavelength 
observations, such as the Spitzer 160 //m band, are observed 
to be correlated with each other ( |Bendo et al.|2008] l, and are 
thought to be more associated with the cooler diffuse ISM. 
While the diffuse ISM is heated by the radiation from star 
forming regions, it is also heated by the diffuse radiation 
field from older stars, meaning that these bands are not as 
strong SFR traces, especially at low SFRs. 

6.3 Fitting the full UV to FIR SED 

Fitting a self-consistent model over the divide between stel- 
lar and dust emission in the SED is of course one of the 
ultimate goals of SED fitting. However, while modelling ef- 
forts are well-underway (see Section |2.2.3| l, unfortunately, 
few authors have attempted to apply these to large samples 
of galaxies. Inversion techniques are not applicable here, as 
the problem is highly non-linear with many free parameters 
and therefore time consuming. Even for calculation of a li- 
brary of model galaxies, it is challenging to provide models 
that are sufficiently simple, complete and fast to make this a 
practical possibility. 

We here need to bypass instances where modellers test 
their codes on single galaxies (e.g. 



Silva et al. 1998 



Pope scu et al.|2000| |Groves et al.|2008| . This is of course 
a most valuable and necessary step to make sure that the 
model does bear on our understanding of reality and to fur- 
ther our knowledge of the underlying physics. Most of the 
results from these studies have been presented above in Sec- 
tion^ For the potential novice reader of this manuscript we 



nevertheless emphasize at this point that proper filter convo- 
lution and^ 2 fitting cannot be replaced by by-eye passing of 
spectra through photometric data points. 



Using the model by|da Cunha et al. ( 2008 1, in a follow-up 



paper da Cunha et al. ( 2010 1 have demonstrated the strength 
of fitting the full SED from UV to IR wavelengths. By fit- 
ting the full SEDs of ~ 3000 galaxies with GALEX, SDSS, 
2MASS, and IRAS data, they were able to determine the star 
formation rate, the star formation history as measured by the 
specific star formation (t^ s ^current - 10 8 Myr star forma- 
tion rate divided by the past average star formation rate), 
dust and stellar masses and other parameters. They found a 
strong correlation of dust mass (Mj) to star formation rate 
(ip, in M Q yr" 1 ), shown below in equation]^}, as well as re- 
lations between the dust to stellar mass ratio and i^ s , and the 
fraction of IR emission arising from the diffuse ISM and ifr s . 



M d = (1.28 + 0.02) xlOV 



Ul±0.01 



(6) 



This work demonstrates clearly the connection between 
dust mass, star formation history and stellar evolution. 
Iglesias-Paramo ~et al.| ( 2007| have gone to the length of 



using the GRASIL code ( |Silva et ai]|1998| l to compute a 
library of 5000 model galaxies and then use Bayesian infer- 
ence to derive the properties of their sample. Their general 
results agreed well with independent studies by other au- 
thors, thus lending support to the notion that full SED fitting 
is a reliable tool to derive galaxy properties. More impor- 
tantly in the present context, they show that their reduced x 1 
distribution has a median value of 2.6, albeit with a long tail 
extending well above 10. Thus the GRASIL library is found 
to reproduce their sample fairly well. Nevertheless, from the 
point of view of reliable SED fitting tools, a more thorough 
analysis of the outliers (model uncertainties, incomplete li- 
braries, AGN, etc.) would be valuable, not only in this but 
in many other works. 



Noll et al. ( 2009 1 present a new code, which they call 



"CIGALE", which effectively computes a library of model 
galaxies and then uses a modified version of the Bayesian 



inference described in Section 4.5 to determine the galaxy 
properties. Diagnostic plots like their Figure 14, which 
shows the residuals between best fit model and data for their 
full sample, are a very useful tool to understand model sys- 
tematics. In their case for example, they conclude that "For 
MIPS 160 jjm the significant deviations can partly be ex- 
plained by the lacking flexibility of the one-parameter mod- 
els of |Dale and HeIoul ( ^2002l l. , ' 

As shown by the last two examples, most UV-FIR SED 
fitting codes are still in their testing phase and have mostly 
been used to confirm results already obtained from more tra- 
ditional single-tracer analyses. The large number of derived 
parameters and our still limited knowledge of their respec- 
tive degeneracies and systematic uncertainties make it dif- 
ficult to go a step further and fully use the full power of 
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SED fitting. Indeed, for the moment it is still question- 
able whether it is not more fruitful to use a combination of 
single tracers to derive one property well (e.g. |Ken nicutt 
|et al.|2009| for SFR). On the long run, however, SED fitting 
holds the promise to provide a large set of galaxy proper- 
ties for large samples. Self-consistent inter-comparison of 
sub-samples with different properties, such as masses and 
SFRs, and the exploitation of constraints on hitherto uncon- 
strainable parameters, such as the relative weights of young, 
intermediate age and old populations, are an exciting avenue 
to explore further in the future. 

7 Conclusions 

We have presented an overview of some of the achievements 
and challenges related to fitting the Spectral Energy Distri- 
butions of galaxies. SED fitting can be used effectively to 
derive a range of physical properties of galaxies, such as red- 
shift, stellar masses, star formation rates, dust masses, and 
metallicities, with care taken not to over-interpret the avail- 
able data. To allow for more progress in galaxy evolution 
studies from SED fitting, we suggest two main areas. On 
the one hand there still exist many specific issues such as es- 
timating the age of the oldest stars in a galaxy, finer details of 
dust properties and dust-star geometry, and the influences of 
poorly understood, luminous stellar types and phases. The 
challenge for the coming years will be to improve both the 
models and the observational data sets to resolve these un- 
certainties. On the other hand, the robustness and accuracy 
of SED-fitting-derived properties still need to be assessed 
more completely. The challenge here is to develop and un- 
derstand the interplay between the fitting routines and the 
available data and models. 

In the hope of accompanying these challenges, the 
present review will be made available on a webpage (sedfit- 
ting.org) together with links to relevant models, fitting codes 
and datasets. We would like to encourage the community to 
send in suggestions for additions and changes to the texj^] 
through this webpage. The intention is twofold: 1) We hope 
to bolster the information currently available in this review 
and keep it up to date over the coming years. 2) Due to our 
bias to the workshop participants we did mention many im- 
portant works in passing, or indeed missed them. We there- 
fore hope that particularly those members of the community 
whom we missed will take to opportunity to add their part of 
the story, thus expanding the current text beyond its original 
scope. 
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Fig. 17 Rest-frame residuals (Am = {template — 
galaxy) I template) for all galaxies with reliable redshifts 
from the zCOSMOS sample. The feature with the strongest 
significance is at around 3500 A, where the templates are too faint 
by 0.05 mag compared to the data. Additional nebular emission 
was added to the original SPS templates for this plot. [Courtesy P. 
Oesch] 
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Fig. 18 IRAC [3.6]-[4.5] color-redshift plot for a sample of 
GOODS galaxies. The solid curves show the change in colours 
with redshift of five different empirical template SEDs. [Courtesy 
M. Brodwin] 
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Fig. 21 Cumulative fraction of objects with 6 Z smaller than agiven 
value. Red line indicates the simulation in which all galaxies have 
been set to have SN=250 in all BVRz; orange indicates a simula- 
tion ith SN= 1 00; and so on. Right panel shows all galaxies, and left 
panel shows galaxies with ODDS > 0.9. Note that only 6.4%and 
1.2% respectively of objects with SN=10,5 have ODDS > 0.9 
(figure from Margoniner and Wittman 2008 1 [Courtesy V. Mar- 
goniner]. 



Fig. 19 Zf-band residuals vs. color and redshift in a sample of 
GOODS galaxies. An effective offset of ss 0.3 mags is apparent. 
[Courtesy M. Brodwin] 
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Fig. 20 Iterative improvement in photometric redshift estimation 
via this simple calibration technique. [Courtesy M. Brodwin] 




Fig. 22 The probability density as a function of the redshift for 
early- and late-type galaxies {upper and lower panels, respectively) 
at different distances marked by the vertical dashed lines. For every 
object, the dotted line shows the empirical relation of p{z\x = m q ), 
and the solid line illustrates the final result of p(z\y = m q , M) after 
properly folding in the photometric uncertainties via the mapping 
in the model (Budavari 2009a figure from) [Courtesy T. Budavari], 
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Fig. 23 Monte Carlo library of 50,000 models created from the 
Version 2007 of |Bruzual and Charlot| ( [2003) . The (g-i,i-H) space 
is binned in cells of 0.05 mag x 0.05 mag and marginalized over 
M/L in each cell. The right panel shows M/L H as a function of 
(g—i,i—H). The scatter at each position is typically 0.11 dex, with 
peaks at ~0.2 dex in the bluest corner. The left panel shows M/L 
as a function of g—i (median value marginalized in a given g—i 
bin). The minimum and maximum M/L values derived from the 
right panel at given g—i are displayed by dashed lines. Figure from 
|Zibetti et alT] < |2009| > [Courtesy S. Zibetti]. 




Fig. 24 Panel a) shows the color-composite RGB image (H in red 
channel, i in green and g in blue) after adaptive smoothing. Panel 
b) plots the distribution of pixels in color-color space (log inten- 
sity scale). The cross marks the position of the global colors of the 
galaxy. The over-plotted contours show the number density distri- 
bution of models from the library described above. Interestingly, 
observed colors are reproduced by the models over the full range. 
Panel c) shows the mass-to-light ratio map derived with the method 
described above. Finally, Panel d) displays the stellar mass surface 
density map, obtained by multiplying the H band surface bright- 
ness times M/L. Images a) and d) show that the galactic structure 
is much smoother in the mass image, where spiral arms have much 
less contrast with respect to what the light shows. Figure from |Zi-| 
|betti et al.|j2009] > [Courtesy S. Zibetti]. 




Wavelength (jim) 

Fig. 25 An example of combining Spitzer-IRS data with a PAH- 
FIT mid-infrared model. The IRS data arise from the central region 
of M82 and are marked by the boxes. The PAHFIT final spec- 
trum passes through these points (green line) and can be seperated 
into PAH features (broad blue curves), fine-structure emission lines 
(purple narrow cuves), hot dust continuum (underlying red cuves) 
all experiencing broad dust extinction, especially the pronounced 
silicate 10 and 18 /im features (represented by the dotted curve). 
Data are from ( Beirao et al. 2008 ). 
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Fig. 26 The dependence of galaxy star formation history as a 
function of distance from the starburst relation of Calzetti ( 19971; 
|Meurer et al.] l [l999) . The lefthand axis is an observable diagnostic 
of the birthrate parameter, the current star formation rate normal- 
ized to the average star formation rate. The righthand axis shows 
the number of years (continuous) star formation has been occur- 
ring, as measured from theoretical spectra. 
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